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Abstract 

A simple pedagogical introduction to the effective action method of quantum 
field theory is given at a level suitable for beginning postgraduate students. 
It is shown how to obtain the effective potential at zero temperature from a 
regularized zero-point energy. The results are applicable to curved as well as to 
flat space. The generalization to finite temperature is also given. It is shown how 
to obtain high temperature expansions of the thermodynamic potential for the 
neutral free Bose gas and the charged Bose gas in both the relativistic and non- 
relativistic limits. The results are obtained for an arbitrary spatial dimension 
and in curved space. Results arc also obtained for the self-interacting relativistic 
gas in three spatial dimensions. A detailed discussion of how the formalism may 
be applied to study Bose-Einstein condensation is given. The interpretation of 
Bose-Einstein condensation as symmetry breaking is discussed. Application is 
also given to the study of charged bosons, both relativistic and non-relativistic, in 
a constant magnetic field. The Meissner effect is obtained for the non-relativistic 
Bose gas in three spatial dimensions. The final application is to the study of non- 
relativistic bosons in a harmonic oscillator potential trap. A number of analytical 
approaches to this are discussed. 

Foreword 

These notes represent the lectures given at The 15th Symposium on Theoretical 
Physics August 22-27, 1996 at Seoul National University. The intent of these lec- 
tures was to present an introduction to the use of the effective action method in 
quantum field theory at a level appropriate to a beginning postgraduate student. I 
assume that the prospective reader has at least a nodding acquaintance with canonical 
quantization in quantum field theory but no detailed knowledge of the more advanced 
aspects. I purposely avoid the use of path integral methods and the formal properties 
of the effective action since they are not necessary for the examples I will discuss. 
It is hoped that these notes will make it easier for a student to make the bridge be- 
tween the elementary concepts of quantum field theory and the more advanced ideas 
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used in contemporary research. Particular emphasis is placed in these lectures on the 
application to Bose-Einstein condensation (BEC), a topic of much current investiga- 
tion. A much more comprehensive introduction to the effective action method will 
be contained in a forthcoming book [|TJ . 



1 Action Principle in Classical Physics 

It is well known how the action principle can be used in classical mechanics to de- 
rive the equations of motion and to discuss symmetries. (See Ref. |2| for example.) 
The Euler-Lagrange equations can be found by demanding that the classical action 
functional 

S[x(t)} = 1* dtL(t) , (1.1) 



be stationary under all possible variations of the path x(t). In Eq. |1.1| L(t) is the 
Lagrangian of the system, which is defined to be the difference between the kinetic 
and potential energies. The square brackets in Eq. |1.1| are used to denote a functional, 
which may be regarded as a real or complex valued function whose domain is a space 
of functions. The Euler-Lagrange equations follow from the demand that S[x(t) + 
Sx(t)} = S[x(t)} to first order in the infinitesimal variation of the path 8x(t). If instead 
of the Lagrange form of the equations we are interested in the Hamiltonian form, then 
we may take 

S[x(t),p(t)} = f * dt{p(t)x(t) - H(t)} , (1.2) 

and demand that it be stationary under independent variations of the path x(t) and 
the momentum p(t). H(t) is the Hamiltonian for the system. (See || for more 
details.) 

The action principle extends immediately from classical particle mechanics to 
classical field theory. The prototype is the vibrating string in one spatial dimension. 
Let ip(t, x) represent the displacement of a point on the string at position x at time t. 
If we consider a small section of string of length d£ and assume that the displacement 
of the string from its equilibrium position along the x-axis is small then d£ ~ dx. This 



I , s f dip(t, x) 
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little section of string has kinetic energy -(pdx) I — — - — J and potential energy 
1 / dip(t,x)\ 2 

-(rdx) I — — - — ) where p is the mass per unit length and r is the tension. The 



2 V ' \ dx 
Lagrangian for this system is therefore 




2 



where the integration extends over the total length of the string. The integrand of 
Eq. |1.3| is referred to as the Lagrangian density. Variation of the action defined as 
in Eq. |1.1| with L(t) given in Eq. [L3] with respect to an arbitrary variation 5ip(t,x) 
leads to the familiar one-dimensional wave equation. The extension to any number 
of spatial dimensions is obvious. 

As already mentioned, one advantage of using an action principle approach is that 
it is possible to discuss symmetries of the theory in a simple way. This is done by 
the requirement that the action functional be a scalar under the symmetry transfor- 
mation. For example, for a particle moving on a curved background (eg. a particle 
constrained to move on a surface) the action should be a scalar under a general 
coordinate transformation on the curved background. This reflects the physical re- 
quirement that the motion of the particle should not depend on the arbitrary choice 
of coordinates used to describe the curved background. If we consider relativistic 
field theory, then the action should be a scalar under Lorentz transformations for 
a theory in Minkowski spacetime. For fields in curved spacetime, the action should 
be a scalar under local Lorentz transformations as well as under general coordinate 
transformations. For gauge theories we would demand that the action be a scalar un- 
der the action of the gauge group. The requirement of invariance leads via Noether's 
theorem to a conservation law. (See || for example.) 

2 Action Principle in Quantum Physics 

Because classical physics is assumed to be an approximation to the more fundamental 
quantum physics, it is natural to ask if there is a quantum analogue to the classical 
action principle. In particular we can ask the question Is there a quantum version of 
the classical action functional ?. One approach to this was initiated by Schwinger in 
the early 1950's in a series of papers ||]. In this work it was shown how a quantum 
version of the classical action principle could be formulated, and which led to the 
Heisenberg equations of motion for the quantum operators. Instead of pursuing this 
formalism directly, we will proceed in a more heuristic manner. Put simply, we wish 
to find an action functional T[<p] for the basic field variables (p which includes all 
possible quantum corrections to the classical action S[(p]. Thus we must have 

r[<p) = s[<p) + o(h), (2.i) 

where the terms of order h denote the quantum corrections to the classical theory It 
is customary to call T[(p] the effective action. 

The effective action would be expected to share as many symmetries of the classical 
theory as possible, although one should be aware that some classical symmetries (such 
as chiral or conformal invariance) may be broken by quantum effects. (The physical 
reason for the breaking of some classical symmetries in the quantum theory is that in 
quantum field theory the process of regularization necessarily introduces a mass scale 



3 



into the theory; therefore any classical symmetries which depend on the property of 
masslessness in the classical theory may be broken.) Particular symmetries which 
we would expect T[<p] to enjoy are Lorentz invariance if the theory is in Minkowski 
spacetime, or general coordinate and local Lorentz invariance if the theory is in curved 
spacetime. In general we would also expect F[(p] to be gauge invariant if the classical 
theory is gauge invariant. 

Although it is possible to be more general, we will limit our attention to spacetimes 
which are ultrastatic ||. In this case the line element may be written as 

ds 2 = dt 2 - g lj (x)dx i dx j , (2.2) 

with t the time and x % {% — 1, . . . , D) the spatial coordinates. We will assume that 
the spacetime Ai ~ R x E with E a D-dimensional Riemannian manifold which is 
compact with or without boundary. In the case where we are interested in Minkowski 
spacetime, rather than dealing with a space of infinite volume, it often proves con- 
venient to restrict attention to a finite box with some boundary conditions, usually 
periodic, imposed on the fields on the box. Our assumption that M. ~ R x E is just 
a generalization of this procedure. In general if we are interested in the case where 
E is really non-compact, we will assume that we have added some spatial boundary 
to obtain a compact spacetime with boundary, and then look at what happens when 
the boundary is removed to infinity at the end of any calculation. The restriction to 
static spacetimes is to ensure that there is an unambiguous definition of a vacuum 
state, and what is meant by the notion of a particle. 

Given the line element in Eq. [2.2| the invariant volume element for the spacetime 

is 

dv x = a/ g(x) dtd D x , (2.3) 

where g(x) = det gij (x) . We will often use 

At x = v^M d D x , (2.4) 

to denote the invariant volume element on E. 

The classical action is defined in terms of the Lagrangian as in Eq. In some 
cases we may wish to define the Lagrangian density £(t, x) in terms of the Lagrangian 
by 



L(t)= dcr x C(t,x), (2.5) 

so that 

%>]= [ dt fda x £{t,x) . (2.6) 
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For the effective action V[tp] we can define an effective Lagrangian density in a com- 
pletely analogous way : 

r>] = [' dt fdcr x C cS (t,x) . (2.7) 
Jh Jt. 

Given an expression for T[ip\ we can obtain the effective field equations from the 
requirement that Y\ip + 5ip] = T[tp] to first order in the arbitrary variation Sip. 

Because the Lagrangian is the difference between the kinetic and potential en- 
ergies, we can define a kinetic energy density and a potential energy density in an 
obvious way. In many situations of interest the potential energy contains the interest- 
ing physics. This is the case if the field is assumed to be constant, or more generally 
independent of time. In such a situation we can write 

r[<p] = -fo-t 1 )VvV eS , (2.8) 

where Vs is the volume of S, and V e g is the effective potential. (V e s is really the 
effective potential density, but it is customary to just call it the effective potential.) 



3 Massive Real Scalar Field 

The classical action functional describing a massive real scalar field is 

S[<p] = jT 2 dtjda x {-3"^ - ~mV - ~£% 2 - £%)} ■ (3.1) 

Here R is the scalar curvature of E, m 2 is the mass associated with the field, U(ip) 
represents the potential for any self-interaction, and £ is a dimensionless coupling 
constant which couples the scalar field to the curvature. Setting the variation of S[ip] 
to zero results in 

Dtp + m 2 (p + £Rfp + U'(tp) = , (3.2) 

d 2 1 
where □ = — — — V 2 with V 2 = —pd^^fgg^dj) the covariant Laplacian on S. 
ot g 

We will initially assume that U'(ip) = 0, so that U(<p) is an arbitrary constant 

which does not affect the field equations. It is impossible to solve the classical field 

equation |3.2| exactly. However we may proceed in a way which will prove useful later. 

Assume that we have the solutions to 

(-V 2 + £i?)/ n (x) = a n /„(x), (3.3) 

and that the set of solutions {/ n (x)} is a complete set. The index n on / n (x) and the 
eigenvalue a n is a multi-index which labels the different solutions. We can normalize 
the solutions by 

^da x / n (x)/*,(x) = 5 fin /. (3.4) 
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It is easy to see that a complete set of solutions to Eq. is then obtained from 



y; n (t,x) = e-^7n(x), (3.5) 

where 



w n = V °n + m? . (3.6) 

The classical contribution to the effective potential is obtained very simply from 
Eq. |3J] with U(ip) = c, where c is a constant, to be 

K ( ff 0) = ~(™ 2 + £R)<P 2 + c (3.7) 

If we wish to calculate the quantum corrections to this classical result we merely need 
to compute the energy density associated with the energy levels hu n with oj n given in 
Eq. j3.6[ Because the scalar field may be considered as a collection of simple harmonic 
oscillators labelled by the multi-index n, we have 



since each oscillator has a zero-point energy of \hu) n . Because uj n will not be bounded 
as n — > oo, the sum over zero-point energies in Eq. |3.<j| will diverge, and therefore 
must be regulated. One beautiful way to do this is by (- function regularization || 0] . 

The basic idea behind (-function regularization is to define the divergent sum 
in Eq. |3.8| by analytic continuation of a convergent sum. We may define an energy 
^-function cu(s) by 



UJ(S) 



J2^n{^n)~ S , (3.9) 



where s is a complex variable, and £ is a constant with units of length introduced 
so that the energy (-function uj(s) has the same dimensions for arbitrary s as the 
original sum in Eq. |3.8| . The idea now is to use the unboundedness of u n as n — > oo 
to our advantage by evaluating uj(s) in a region of the complex s-plane where the 
sum converges, and then evaluating co>(0) by analytic continuation of the result. We 
therefore define the regulated expression for \Qg by 

= ^(0) . (3.10) 

It may be that uj(s) has a pole at s = 0. In this case uj(0) denotes the analytic 
continuation of uj(s) back to a small neighbourhood of s = 0. We can give a very 
loose justification for assuming that there is a region of the complex s-plane where the 
energy (-function converges as follows. (We will see this in a more general way later.) 
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If we consider the behaviour of u n for large n, then we can approximate the discrete 
set of eigenvalues a n with the continuous value of k 2 . Then the convergence of u>(s) is 
determined by the convergence of / d D k(k 2 + m 2 )^' 2 . This integral is convergent 
for 5J(s) > D + 1 and we therefore expect that the energy (^-function Eq. |3.9| converges 
in this region of the complex plane. 



3.1 Free Scalar Field in Flat Space 

It is helpful to illustrate how this regularization works in some simple examples. We 
will first consider the scalar field in a box in flat spacetime with periodic boundary 
conditions imposed on the box walls. If L±, . . . ,L D represent the sides of the box, 
then the eigenvalues a n defined in Eq. EO are just 



' 2-nn^ 



i=i 



Li 



(3.11) 



with rii = 0, ±1, ±2, . . . and n standing for the D-tuple (ni, . . . ,tid). The energy 
^-function in this case is 



UJ(S 



£•■■£' 

ni=— oo ™<i=— oo 



£ 



i=i 



+ m 



(1-.J/2 



(3.12) 



If we are interested in the case where the sides of the box become very large, ultimately 
we will take the infinite volume limit, then we can approximate oj(s) by replacing the 
sums over rii with integrals. (We will return to this at the end of the lectures.) 
Changing variables of integration leads in a straightforward way to 



vis) 



Vv 



d D k 



(2tt) d 

The integral may be evaluated using the identity 

1 



k 2 +m 2N(l-.)/2 



(3.13) 



r z) 



z — 1 —at 



dt t z ~ L e 



(3.14) 



which is valid for 3l(z) > and 3?(a) > 0. This gives 



1 



d D k 



T((s-1)/2)J (2^J 
1 



du (s-3)/2 e -(k2W)t 



(4tt)- d / 2 / dt&- 3 - D ^ 2 e- m 1 



r((s-i)/2) 

^ r((g l _ 1)/2) (47r)^/ 2 r((, - 1 - D)/2){m 2 r+^l 2 . (3.15) 
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The middle line follows from performing the integration over k, which just involves a 
product of gaussian integrals, and the last line follows from using the identity Eq. |3. 14 



again. We must assume that dt(s) > D + 1 to justify the steps leading up to Eq. |3.15| . 

The result for our energy ^-function in Eq. |3.15| can be seen to have simple poles 
whenever (s — 1 — D)/2 — —j for j — 0, 1, 2, . . . coming from the poles of T((s — 1 — 
D)/2). We can analytically continue the result in Eq. p. 15 throughout the complex 



plane. a>(s) is analytic at s = if D is even and has a simple pole at s = if D is 
odd. We can examine the cases of even and odd dimensional spaces separately. 
For D even we have 

.(0) = V- s (4.)^ r( - ( ( 1 ^ / ^ /2) (^)^)/ 2 , 



which leads to 



as the quantum correction to the classical potential. (The result has been simplified 
using some properties of the T-function.) 

For odd D we must expand about the pole of the T-function. It may be shown 
that 

V '« ~ 2 (fi±l)! \ in) Is + 2 + 3 + 

is the regularized expression for the quantum correction to the classical potential. 

There is one obvious difference between the results found for even and odd dimen- 
sions. When D is even, the regularized expression for is finite at s — 0, whereas 
for odd D the result diverges as s — * 0. In order to deal with the divergent term we 
must renormalize the effective potential. This can be done very simply in the present 
case by imposing the renormalization condition 

Kff(y = 0) = . (3.18) 

This leads to 

= V cS ^ = 0) = V^^ = 0) + V^(cf = 0) 

= C + Vg> (3.19) 

since does not depend on ip in the present case. The renormalization condition 
results in the constant term in the effective potential cancelling the quantum part of 
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the effective potential. For D odd, this also removes the divergent part of Eq. |3 . 1 7| . 
We are left with the final expression 

V eS {cp) = ^mV , (3-20) 

as the renormalized effective potential. In this case the result is not particularly 
interesting. 



3.2 Casimir Effect 

A slightly less trivial example occurs if we consider a flat space with some boundary 
conditions imposed on the field. A prototype for such an example is the Casimir effect 
H which is concerned with vacuum fluctuations in the electromagnetic field between 
a pair of neutral parallel conducting plates. As a scalar field analogue of this we will 
consider the massive scalar field which satisfies Dirichlet boundary conditions on two 
parallel planes at x 1 = and x 1 = L\, but which is unconfined in the remaining spatial 
directions. We will choose the spatial dimension to be D = 3 here for simplicity of 
presentation. The mode functions which obey Eq. |3.3| (with R = again) may be 
chosen to be 



(3.21) 



3=1 



where A n is a normalization constant, n\ = 1, 2, ... , and rij = 0, ±1, . . . for j = 
2,3. We will look at the limit L 2 ,Ls —> 00, but will keep Li, which represents the 
separation of the plates, finite. Periodic boundary conditions have been imposed in 
the directions transverse to the plates. Note that we cannot allow n\ — because 
this value does not lead to a normalizable eigenfunction. Negative values of n\ can 
be neglected because the sign of the normalization constant is free. 
We now have 



CTr, 



nrii \ 2 ( 27rn, x 

j=2 V J 



Li 



(3.22) 



in place of Eq. |3.11| . From Eq. |3.9| we have 



0J(S 



ni=l ri2=— 00 n3=— 00 



2 3 /o x 2 
7T711 \ v"^ / 27m,- 

3=2 



+ m 



(l-«)/2 



If we take the limit L 2 ,L 3 — ► 00, and replace the sums over n 2 and n 3 with integrals 
as before, it is easy to show that 

(3-s)/2 



wis) 



L 2 L 3 T((s - 3)/2) 
4tt r((s-l)/2) 



E 

ni=l 



m 



rrrii 



(3.23) 
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It is not possible to replace the remaining sum over n\ with an integral as we did for 
the sums over n 2 and n% because we are not interested in the limit L\ — > oo. There 
are a number of ways to evaluate the sum appearing in Eq. |3.23| in terms of a contour 
integral ||. Instead of trying to obtain an exact result, we will look at the expression 
in the limit mL\ « 1. This limit is of interest if we wish to study the massless limit. 

When mLi « 1 we can expand the summand of Eq. 3.23| in powers of mLi in 
an obvious way. If we keep only the first three terms, we have 

3~s oo 



.LoL 



UJ(S 



2-^3 



4tt 



1 



1 

+ 2 



ni=l 

mLi 



mLi 



71 



717 



7i 



7l-\ 



LoL 



2-^3 



47T V L 1 



7T 



3-s 



[8-3) 



Cr(s - 3) 



mLi 



7[ 



1 

5' 



1-s 



mLi 



71 



Cr(s + 1) 



(3.24) 



The definition of the Riemann ^-function (r(s) = s for 3?(s) > 1 has been used 



n=l 



here. The only pole as s — > arises from the term involving Cr( s + 1) which has a 
simple pole. Using properties of the Riemann (^-function fH| we have (noting that 
Ve = LiL 2 L 3 ) 



v; 



(i) 



off 



hTT 2 

8Lf 



180 



1 

12 



mL\ 



71 



mL\ 



71 



4 r 



+ 7 - X^tjLx 



(3.25) 



If this result is compared with that in Eq. p. 17 it can be observed that the pole parts 
of both expressions are the same. Thus the divergent part of the effective potential 
is unaffected by the boundary conditions imposed on the field. This is a special case 
of a more general result JTTJ] which we will establish in the next section. A natural 
thing to do here is to choose the constant c to be the same as that found for the 
same theory in Minkowski spacetime. From Eqs. [3.17 and 3.19 we find that the 
renormalized effective potential is 

1 



2 2 

m (p 



hm 4 
'32^ 



n 2 h 

1440L? 

1 1 , 

- + 7 + - In 

4 [ 2 



15 



mLi 

71 

m 2 L\ 

ATI 2 



In the massless limit we are left simply with 

7T 2 h 



off 



1440L? 



(3.26) 



(3.27) 



10 



which is half of the result found for the electromagnetic field. 



3.3 Xcp 4 Theory 
Suppose that we take 

A 



[%) = i y + c. (3.28) 



The classical field equation |3]2] becomes 

Dtp + m 2 ip + £Rip + -cp 3 = . (3.29) 
6 

In order to use the energy ^-function method we must know the energy levels as- 
sociated with quantum fluctuations about the classical solution to the equation of 
motion. To this end, write 

<p = tp + iP, (3.30) 

where if) is the quantum fluctuation in the field. If we are only interested in work- 
ing out the lowest order correction to the classical effective potential, then we may 
substitute Eq. |3.30| into Eq. [3.29| and linearize in ip. This results in 

Dip+m 2 ip + £Rip + ^ 2 ip = . (3.31) 

We have not assumed that <p is constant here, but in conformity with our assumption 
on the static nature of the spacetime we will assume that (p = y?(x) only. 
The classical potential is 

V£\<p) = c+ l -\V^ + \m 2 p 2 + l -HRf + (3 32) 

In place of Eq. K73], let 



-V 2 + iR + ^<p 2 ) /„(x) = <r n / n (x) . (3.33) 

Take ip n (t, x) = e~ luJnt f n (x) as in Eq. |3.5| , so that uj n is again given by Eq. p.6[ . We can 
still define V$ by Eq. |XTD| with uj(s) given by Eq. p.9| . For general S and arbitrary 
ip it is impossible to solve Eq. |3.33| explicitly for the eigenvalues a n ; nevertheless, we 
can apply a very general method to obtain approximate results. 

First of all we will examine the divergent part of the effective potential in curved 
space to study the renormalizability of the theory. From Eq. |3j| we can write 
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Using Eq. |3.14| we can obtain 



where we have defined 

0(t) = J2e~ t{an+m2) . (3.35) 

n 

The advantage of this procedure is that a great deal of information is known about 
@(t). Our first observation is that 0(t) t — > oo where a is the 

smallest eigenvalue. Assuming that <7n > this shows that the integrand of Eq. |3.34 



is well-behaved as t — > oo for all s. This means that if u}(s) has any poles, they can 
come only from the t — > limit of the integrand. We may split up the integration 
range in Eq. |3.34 into two parts as follows : 



uj(s) = ujp(s) + u r (s) , (3.36) 



u, (s) = — - [ t0 dt&- 3)/2 &(t) , (3.37) 



with 



T((s - l)/2) 7 

■* ( " = f((73T W i i /" , *- 3,/2 eW. (3-38) 

From our above discussion because of the way we have split up the integral, u r (s) 
should be an analytic function of s, and any poles in u(s) should come only from 
u)p(s). Here t is an arbitrary positive constant which we can specify freely. It may be 
taken to be arbitrarily small, so that to see if u>(s) has a pole at s = we only need 
to know the behaviour of 0(i) for small t. Fortunately this is a well studied problem 
I2| . It can be shown that 0(i) has the asymptotic expansion 



9(t) ~ (Ant)- D / 2 Yl tk0 k , ( 3 - 39 ) 
fe=o i,i,„. 



for some expansion coefficients 9k which are independent of t. Using this in Eq. |3.37 
results in 



,fc+0-D-l)/2 
En 



-(47r)-^ 2 V 9 k —^ — 7T77T i (3.40) 



r((s-l)/2) v ' "k + (s-D-l)/2 

' 2 ' ' 



where we initially assume dt(s) > D + 1. (This establishes our earlier claim that ou(s) 
is analytic in s if 3?(s) > D + 1 as stated earlier.) Analytic continuation back to s = 
shows that uj(s) has a pole only if k — (D + l)/2. We have 

cu(s) = -(4nr {D+1)/2 9(D + i)/2 - s + ■ ■ ■ , (3.41) 
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in a neighbourhood of s = 0. Keeping only the pole term, we find 

= -^r (D+1)/2 d(D +m (3.42) 

All that we require is a knowledge of the coefficients 9k appearing in Eq. [3.39| . 

For any operator of the form (—V 2 + Q), where Q is any arbitrary function of x, 
the first few coefficients are known for £ any smooth compact Riemannian manifold 
with smooth boundary. If we concentrate on manifolds with no boundary, then 

(3.43) 

l R -0) > ( 3 - 44 ) 



9 = 




0i = 






h 


9 2 = 






h 



1 1 „o 1 - 1 



— V 2 i? + —R 2 — ——RijR 1 -' + ——RijkiR % ^ kl 



30 72 180 J 180 



--^RQ + ^Q -q^Q) ■ ( 3 - 45 ) 

The expressions for 9 3 and # 4 are known but are too lengthy to give here. If £ has 
no boundary, then #1/2, #3/2, • • • all vanish. From Eq. |3.42j this means that there is no 
pole term when D is even. For the physically interesting case of D = 3 we find 

VcS - 32n 2 V* s + 1 J 

It is straightforward to show that when U (ip) = Xip A / 4! the theory is renormalizable 



in curved spacetime with counterterms in agreement with those of Ref. fllT 



Although we have examined the problem of computing the infinite part of the 
effective potential, for general <^(x) it is difficult to evaluate the finite part. As an 
example we can consider the case where £ = M D with <p constant. In this case it is 
evident from Eq. |3.33| that the presence of (p acts just like an additional term in the 
mass. We may simply obtain from our previous result for the non-interacting 



r(l) 
'eff 

theory in Sec. |3J] by replacing m 2 with m 2 + \<p 2 ■ If we take D = 3 in Eq. |3.17| we 
find 



u 2 + Vu--- 

cS ~ 64tt 2 V + 2M U 2 



2,^-2 

4 + 2^ 



(3.47) 



when we replace m with m + £ip . The total effective potential is 

^^c+imV + ^ + yW. (3.48) 

In order to renormalize the theory we must impose some renormalization condi- 
tions on V e ff as we did for the non-interacting theory. We will write 

c = CR + h5c, (3.49) 

m 2 = m 2 R + hSm 2 , (3.50) 

A = Xr + HSX, (3.51) 
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where the subscript "R" denotes a renormalized quantity. Because is already of 
order h we may take all quantities which occur in it as renormalized. We will choose 
the renormalization conditions 



V cS ((p 



dp 2 
dp> 4 



off 



cff 



0) 



y?=o 



(p=o 



CR 



m 



R > 



A 



R 



(3.52) 
(3.53) 

(3.54) 



The values on the RHS of these conditions are chosen to agree with what is found 
for the classical potential. It is now straightforward to obtain the counterterms and 
hence the renormalized effective potential. The result is 



V eS ((p) 



1 2 -2 , -Vr -4 

CR+^ m RV + 4|-^ 



512tt 2 ^ 



h 



64vr 2 



1287T 

2 , -2 



In 1 + 



2m| 



(3.55) 



If we try to let m 2 R —>■ in Eq. [3.55| to obtain the result for the massless theory 
there is seen to be a problem. This is due to our choice of the renormalization 
condition Eq. |3.54| . Although the presence of <p alters the mass from its tree level 
value of zero to ^p 2 when we impose the renormalization condition at p = there 
can be problems due to infrared divergences. To deal with this we can replace the 
condition Eq. |3.54| with 

= Xr, (3.56) 

(f=M 



dip' 



-Vc 



cff 



where M is an arbitrary renormalization point [|T^]. This leads to the result 

A _„ h\l _ A f, ff_\ 25 

M 2 



V eS {<p) 



-4 fc^R 



In 



6 



(3.57) 



This expression is usually referred to as the Coleman- Weinberg effective potential. 
Ideally we would now solve the effective field equation V^ s (<p) = for (p. It is obvious 
that <p = is one solution to the effective field equation. Unfortunately as discussed 
in Ref. [13| other solutions lie outside the range where perturbation theory is valid, so 



that there are no non-trivial self-consistent solutions to the effective field equations 
apart from ip = 0. 

The fact that V e g(<p) in Eq. |3.57| depends on an arbitrary renormalization constant 
M might seem to be a weakness in the result; however it is actually a great strength 
since it leads to a study of the renormalization group [13]. If a different choice is made 



for M then it is necessary to redefine the renormalized parameters in the theory in 
a way which makes all physical predictions unchanged. We will leave this interesting 
subject as outside the main aim of the lectures. 
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4 Some Thermodynamics and Statistical Mechan- 
ics 



So far we have ignored the effects of a non-zero temperature on the effective action. 
We now wish to see in a simple way how to include finite temperature corrections to 
the results obtained in previous sections. We will only consider the case of thermal 



equilibrium here. More details can be found in the standard references |L4|, [15], [16 . 

We have evaluated the effective potential at T = by summing up the zero- 
point energies. The simplest way to include the effects of finite temperature is to 
appeal to standard thermodynamics. Suppose that we have a system described by a 
Hamiltonian operator H which has a conserved quantity described by an operator Q. 
Q and H necessarily commute. For a non-relativistic system we may be interested 
in a conserved number of particles for neutral systems, or a total conserved charge. 
For relativistic systems the total particle number cannot be conserved due to pair 
creation, but the total charge should be fixed. All of the relevant thermodynamics is 
encoded into the grand partition function 3 defined by 



tr 



(4.1) 



Here f3 = (kT) -1 with T the temperature, and \x is the chemical potential associ- 
ated with the conserved quantity described by Q. The grand partition function is 
regarded as a function of the thermodynamic variables T, V, \i. For curved space and 
background fields, 3 is also regarded as a function of the metric tensor on £ and the 
background fields which are held fixed in the computation. 

In place of 3 it is often convenient to define the g-potential by [jTS] 



g = hi3- (4.2) 

In the limit where the volume becomes very large, 

q = PPV , (4.3) 

where P is the pressure. The connection between statistical mechanics and thermo- 
dynamics can be made from the first law which reads 

dU = TdS - PdV + fidQ , (4.4) 

where Q is the conserved quantity described by the operator Q, and U is the internal 
energy. (Other terms can be present in Eq. |4.4| if the system is, for example, rotating 
or in an external field.) 

A key thermodynamic quantity is the Helmholtz free energy A defined by 

A=U-TS. (4.5) 
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Computing dA and using the first law in the form of Eq. O results in 

dA = -SdT - PdV + ndQ . 



(4.6) 



This demonstrates that A is regarded as a function of T, V, Q. The importance of A 
lies in the fact that in thermal equilibrium for fixed T, V, Q the Helmholtz free energy 
must be a minimum. (See for example jl7|, [Eg].) 

Finally, the grand potential (sometimes just called the thermodynamic potential) 
ft is defined by p| 



ft = U-TS-nQ 
= A-fiQ . 



Computing dQ and using Eq. 

dn 



gives 

SdT - PdV 



Qdji 



(4.7) 
(4.8) 



(4.9) 



showing that ft depends on T, V, \i. The connection between macroscopic thermo- 
dynamics and microscopic statistical mechanics is made by noting from Eq. [4.7| and 
the non-differential form of Eq. ^4], U = TS — PV + fiQ, which follows from Euler's 
theorem on homogeneous functions and the fact that S, V, Q are extensive variables 



Thus 



ft = -PV . 



q = -/3ft . 



In terms of the energy levels E n of the system we have 



E 

n 



In [1 



From Eq. [4.9| we have 



Q 



l 





dn) 

dq 
d/j, 



T,V 



T,V 



(from Eq. 4.11) 



(4.10) 
(4.11) 
(4.12) 

(4.13) 

(4.14) 
(4.15) 



This is the standard Bose-Einstein distribution function. Another way to show con- 
sistency of all of our definitions is to use the definition of Q as a thermal average 



Q 



-tr 



Qe 



(4.16) 
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and to note from this that 



Q /?3 u!) 



1 fdq 

t,v P\ d ^ 



T,V 



using Eq. $72\ . 

A simple and important example is provided by the simple harmonic oscillator. 
Taking Q to be the number operator we have 



tr (r /3(H-MQ) 



£< 

n=0 

oo 

£■ 

n=0 



n e 



-/3(_E n -/m) 



(4.17) 



where we have defined = E n \n) and = n|n). If we use E n = {n + l/2)u;, 

then the sum in Eq. |4.17| is just a geometric series which is easily summed to give 



3 = e~^ 2 [1 



r 1 • 



Using Eq. |4.11| and Eq. [O] we find 

n 







(4.18) 



(4.19) 



as the grand thermodynamic potential. From Eq. |4.9| we have 

Q = 



_ i] 



-i 



(4.20) 



Note that as T — > (so that /3 — > oo) the second term in Q vanishes and Q 
which is just the zero-point energy. 

In quantum field theory we can consider the theory as an infinite collection of 
simple harmonic oscillators exactly as we did before. If ui n denotes the frequency of 
the mode / n (x) as in Sec. |], then we have 



n 



(4.21) 



The first term appearing on the RHS is observed to take the same form as we had 
in the zero temperature case considered earlier. It is only this part of Q which is 
divergent and requires regularization. We may write 



f2 — Qt=o + ^iyo 



(4.22) 
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where 



n T ^ = - ln I 1 - e^ (Wtt - M) ] , (4.23) 

n 

is the finite temperature correction to the zero temperature theory. The zero tem- 
perature part is 



|l>- (4-24) 



n T=0 2 

n 

From Eq. |3l] this can be recognized as V^g Vs. Thus we can define a finite temperature 
effective potential by 

V£> = £■ ■ (4.25) 

There is a natural split of this into a zero temperature part and a finite temperature 
correction to the zero temperature part. Finally we note that we require \i < uj to 
ensure that flr^o is rea l ; an d to avoid a negative occupation number. 



5 Free Bose Gas in Flat Space 

5.1 Non-relativistic Gas 

We will examine the case of a single real degree of freedom for which the dispersion 
relation is non-relativistic. Here 



J .i 



if we resort to the familiar device of imposing box normalization with the infinite box 
limit taken. If we do not have any conserved quantities then 

f2 T#0 = I / dn x ---dn D ln [l - e~^\ , (5.2) 



after the large box limit is taken. Making a change of integration variable rij — *> 
and expanding 

OO 

ln(l-^) = -Ep ( 5 - 3 ) 

k=l 
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gives 



n 



1 f d n n 



oo 

(3 f^kj (2ir) D 



3 fc^n 2 

e 2m 



after performing the sum over k. If the T = contribution is ignored we have 



(5.4) 



v: 



(i) 



cff 



1 / m \ D/ \ n D, 
U 2"' 



(5.5) 



as the thermal part of the effective potential. 



5.2 Relativistic Gas 



For the relativistic gas we have from Eqs. RT6] and |3.11| that 



D 

E 

j=i 



27rnj 

"l7 



m 



1/2 



(5.6) 



This gives the same expression as in Eq. |5.2| in the large box limit. Changing variables 
of integration as before we have 



n 



Vy, 



7yo 



P J (2tt) 



In 



-,fl(n 2 +m 2 ) 1 / 2 



(5.7) 



It is not possible to evaluate this in such a simple way as in the non-relativistic 
case. Instead we will concentrate on just obtaining the high temperature expansion 
introducing a technique which we will use later. 
Begin by expanding the logarithm using Eq. W. 



Vz ^ 1 f d D n fe/3(n 2 +m2) i/2 

n T ^ = - T ^- k ]— D e 



(5i 



It is still not possible to evaluate the integral in simple terms. For high T, specifically 
(3m << 1 or kT » m, we may obtain a useful expansion by using the Mellin-Barnes 
integral representation for the exponential function : 



I t>c+ioo 

/ daT(a)x~ a . 

2m 



(5.9) 



Here cGl with c > chosen so that the integration contour is to the right of the 
poles of the T-function. The proof of Eq. 0] is simple : just close the contour in the 
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left hand of the complex plane. T(a) has simple poles at a = —n for n = 0, 1, 2, . . . 
with residues (— l) n /n\. The residue theorem gives 



e+ioo 



— : / da T(a)x a = 



2?r« ./,_,oo ~ ri! 



-X 



which may be recognized as the Maclaurin series for e 
Using the representation Eq. |5.9| in Eq. |5.8| leads to 



Noting that 



" 71 ^2 , _2\-a/2 _ ( A „\-D/2 L \ 2 )_ ( m 2\(p-a)/2 



+ m?)-*/* = (4tt)- u ^ _) a ' (m 2 )^/ 2 , (5.11) 



(2vr)^ v 7 v 7 r(a/2) 



which can be established as in Eq. |3.15| assuming 3?(a) > D, if we deform the inte- 
gration contour in Eq. [5.10| to the right so that c > D, we can interchange the orders 
of integration and obtain 

- - ii ^^r? r(a)ir ^ (i+a) ffiS (m2)< " o,/2 

= ^«^r:«a wu „ )W . ( , 12) 

/? VW 2ni J c _ ioo r(a/2) 



This can be simplified by making use of the duplication formula for the T-function 

(5.13) 

We find the following representation for the thermodynamic potential : 

V E fm 2 \ D/2 1 f c+i r ^ fa + l\ ^ (a-D* 



T(2z) = (47r)~ 1 / 2 2 2 T(z)r (z + 



^° = -^UJ 2™ > ' 



c— too 



x( R (a+l)((3m/2)- a . (5.14) 

If we close the contour in the left hand of the complex plane we encounter three 
sources of poles. T (^p) has simple poles at a — — 1 — 2n for n — 0, 1, 2, . . . ; T f^f^) 
has simple poles at a = D — 2n for n = 0, 1, 2, ... ; (r(1 + a) has a simple pole at 
a = 0. Because (r(— 2n) = for ra = 1, 2, . . . it follows that T (^j-) CR( a + 1) has 
simple poles only at a = 0, — 1. It is easy to show that when D is even the integrand 
of Eq. |5.14j has a double pole at a = 0, and when D is odd, there is a double pole at 
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a = — 1. In either case all remaining poles are simple ones. It is therefore expedient 
to examine the cases of even and odd D separately again. 

For even D it is straightforward to show using the residue theorem that 



n 



Vy 



2yo 



OO 



n (D+l)/2gD+l ^ nl 
' n=0 
n^D/2 



Vy / m 



2\D/2 



D + l 
2 

D/2 



n ( R (D + 1 - 2ra) 



-(3 2 m 2 



2(3 V 4tt J (D/2)! 



£ \ - M(3 2 m 2 ) 



k=i 



rrv 



2 x (D+l)/2 
- I 



47T / 



-l-D 



(5.15) 



For odd D 



(D-l)/2 



_ 7r -(D+l)/2 /3 -D-l ^ 



n=0 



nl 



2 \ (D+l)/2 

m \ 
4^) 



p + i 

2 

(D+l)/2 



- n CflP + 1 - 2n) 



(3m 
~2~ 



2n 



fc=l v 7 



E ^(- 2 >— ^|(f) 2 "- ("«) 

n=(D+3)/2 V I ) \ / 

It is possible to simplify the result using properties of the V- and ^-functions. In the 
physically interesting case of D = 3 we find 



x 



n 



Vy 



9(T 24 P 12tt 



m 



64tt 2 



3 rtl //3m 
27 - - + 2 In — 

OO 



V 7 / 2 r 4 V ; , --P 2 m 2 ) . (5.17) 



B n are the Bernoulli numbers. 



6 Charged Bose Gas 

In order to illustrate how the calculations proceed for /i 7^ we will consider the case 
of a conserved electric charge. We will also set up the formalism to deal with a static 
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background magnetic field described by a vector potential A(x). The classical action 
functional for the magnetic field will be written as 

Scm= r dt h a * {\ FijFij ~ jiextA ] ■ (o) 

Here = V^A, — VjAj is the field strength tensor which describes the magnetic 
field. J ext is the externally applied current responsible for setting up the background 
magnetic field. If we want to keep the spatial dimension D general then the magnetic 
field should not be described by a vector field, but rather by the antisymmetric 
tensor field Fij. In the special case of D = 3 we may introduce a vector field B with 
components B l through the duality relation Fij = eijkB k where is the Levi-Civita 
tensor. 



6.1 Non-relativistic Field 

We will consider the non-relativistic theory first. This is described by a complex 
Schrodinger field \I/ and we take the action functional to be 

S[ty, *t] = jf 2 dt Jda* 1 1 (V^ - - 7^H D ^| 2 " ^Wl^l 2 } • ( 6 - 2 ) 

It is easy to verify that under independent variations of \I/ and the Schrodinger 
equation and its complex conjugate are obtained. We have set H — 1 as usual. U (x) 
represents any static potential which might be present. D = V — ieA is the usual 
gauge covariant derivative. We will adopt the gauge choice 

V • A = . (6.3) 

The theory is invariant under the local U(l) gauge transformation 

tf(t,x) -> e iee(x) ^(t,x) , (6.4) 
A(x) -> A(x) + V0(x) , (6.5) 

for arbitrary 0(x). Associated with this gauge invariance is a conserved Noether 
current which can be found by varying the matter part of the action with respect to 
Ai and defining the current by 

5S=- I 2 dt !da^ J' l 5Ai . (6.6) 



The minus sign in Eq. |6.6| is for conformity with that in the current which appears 
in Eq. |6.1| and ensures that the current appears with the correct sign in the Maxwell 
equations. It is found that 

9 

op p~ 

J = — (W^-^W) A|*| 2 . (6.7) 

2m m 
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The charge associated with this current is 

Q = e do x \ty\ 2 . (6.8) 



s 



For this non-relativistic theory we can deal with a fixed particle number rather than 
a fixed total charge. In this case the results of Eq. and ^]8| are divided by e, in 
which case J may be recognized as the probability current of wave mechanics. 
The Hamiltonian for the theory with action Eq. p.2\ is simply 



H = jda x <{ ^|D^| 2 + t/(x)|^| 2 } . (6.9) 



To incorporate the conserved charge we introduce the Lagrange multiplier /i as de- 
scribed in Sec. |] and define 



H = H-fiQ = jda x <J — |D*f - e^|*r + U(x)\V\ 2 } . (6.10) 



Associated with H we can define a Lagrangian L by a Legendre transformation and 
obtain the action functional 

= J dt Jda x | % - (V^ - ^) - ^H D ^| 2 + e/^l^l 2 - ^Wl^l 2 } (6-11) 

which incorporates the fact that the total charge is fixed. We can now use this action 
in the usual way. If ^ represents the background field, assumed to be independent of 
time, it is easily seen that the effective potential density has a classical part which is 
given by 

= + - ^(xPl 2 > (6-12) 

In order to calculate the quantum part of the effective potential we may proceed 
as we have done before and calculate the energy levels for the system. The field 
equation which follow from varying Eq. |6.11| with respect to W is 

= ii> + — D 2 ^ + efjM - C/(x)# . (6.13) 
2m 

Apart from the term in /z which accounts for charge conservation this is just the 
Schrodinger equation. Let {/ n (x)} be a complete set of solutions to 



/„(x) = (T„/„(x) . (6.14) 



Then a n are seen to be the energy levels of the time independent Schrodinger equation. 
If we write 

^(t,x) = ^A n e- li? "7n(x) (6.15) 
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for some coefficients A n , we see that 

E n = a n -efi. (6.16) 
The thermodynamic potential fiiyo is, from Sec. f|, 



n 



^ln[l-e-^] 







^J2 ln [1 - e"^-^] . (6.17) 





This last result could have been written down immediately in a straightforward way 
from Sec. |j; however the same procedure when applied to the relativistic theory is 
less trivial. 

As an illustration, we will evaluate the thermodynamic potential when there is no 
magnetic field present. Then a n is simply the same as u n in Eq. |5.1| . If we ignore e 
in Eq. |6.17 , so that we talk about a fixed particle number rather than a fixed charge, 
then in the large box limit we have 



fiiYo = ^Jdm---dn D In [1 - e-^^] 



(6.18) 



in place of Eq. |5.2| . The presence of the chemical potential complicates things slightly. 
Making use of Eq. Ol we have 



Vr. f d D n 

J In 



oo 1 

J2je~ k ^ n2 -») (6.19) 



The integral over n may be performed as before, but this time the sum over k is not 
simply a Riemann (^-function. We have 

^° = -f(^) (6 ' 20) 

Note that as discussed in Sec. f| we need /x < u = here in order that there be no 
negative occupation numbers. The sum in Eq. |6.20| is therefore convergent. It can be 
written in terms of the polylogarithm function defined as 



oo r 



L ^) = E^> ( 6 - 21 ) 

71=1 

for all p if \z\ < 1, and for 3?(p) > 1 if \z\ = 1. We can write 

D/2 



fi ^o = ~r(£z) Li 1+D/2 . (6.22) 
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6.2 Relativistic Field 

The action for a charged relativistic scalar field $ is 

S = j 2 dt fdcr x |(^$) t (D^) - m 2 ^ - C/(x)$ f $} . (6.23) 

(We will keep Aq temporarily non-zero here.) = — ieA^ is the relativistic 
covariant derivative. The action Eq. |6.23| is invariant under the U(l) gauge symmetry 



$ -> e ied § , (6.24) 
-> A M + <9 M 0. (6.25) 



The Noether current which follows from Eq. |6.23| is 



J M = i e ($t^$ _ 5 M $t$) _ 2ie^$ f $ . (6.26) 
The conserved charge is 

Q = Jda* Jo = ie Jda* (V$ - $ f $) , (6.27) 



after choosing the gauge A = and A satisfying Eq. pT3| . 
The momentum canonically conjugate to $ is 

BC 

n = -Z- = $t (6.28) 
(9$ 

where £ is the integrand of Eq. |6.23| . The Hamiltonian density is 

n = n$ + n f $ f - £ 

= n f n+ |D$| 2 + m 2 |$| 2 + f/(x)|$| 2 . (6.29) 
To incorporate the conserved charge we form 

H = H-fiJo (6.30) 

= n t n+|D$| 2 + m 2 |$| 2 + t/(x)|$| 2 -ie/i(n t $ t -n$) . (6.31) 

H is just H = da^H. 

is 

To go back to a Lagrangian description we form C defined in terms of the Legendre 
transformation of Ti by 

I = n$ + n 1 ^ - n (6.32) 

with 

dl~L 

$ = — = n f + ieu$ . (6.33) 
all 
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Similarly, 



dT~C 

$t = __ = n - ie/i$ t . (6.34) 



If we eliminate II and 11^ Using Eqs. |6.33| and |6.34| we find 

£ = $t$ + i e/ i($t$ _ $t$) + eV* f $ - |D$| 2 - m 2 &$ - U , (6.35) 

as the Lagrangian density which describes the theory with a conserved charge. 
The Euler-Lagrange equations which follow from Eq. |6.35| are 

= -$ + 2ie^ + D 2 $ - (m 2 - eV)$ - C/(x)$ , (6.36) 
= -$t _ 2ze/i $t + D 2 $ f - (m 2 - eV)$ f - [/(x)$ f . (6.37) 

If we let {/ n (x)} be a complete set of solutions to 

[-D 2 + C/(x)] / n (x) = <r n / n (x) , (6.38) 

we may expand 

$(t,x) = ^A n e-^/„(x), (6.39) 



for expansion coefficients A n . Substitution of this expansion for $ into Eq. |6 . 3 6| results 
in 

= E 2 n + 2efiE n -a n - (m 2 - e 2 /i 2 ) , 



giving 



E n = ^a n + m 2 - e/i . (6.40) 

A similar analysis applied to Eq. |6.37| gives 

E n = \fa~ n + m 2 + e/i . (6.41) 

The two results for E n are seen to differ only in the sign of the electric charge. 
Physically they correspond to the contributions from particles and antiparticles. The 
thermodynamic potential in Eqs. [4.22| — 14.24| is therefore 

ftr=o = ^ n , ( 6 - 42 ) 

n 

and 

fi T #o = lj2 ln i 1 ~ e-^-^] + 4 X) ln t 1 ~ e~^ n+etM) ] , (6.43) 

n n 
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where 



\fa n + m? . (6.44) 



Qt=o in Eq. |6.42| differs from that found for the real scalar field by a factor of 2, 



which is easily understood due to the fact that a complex scalar field has two real 
degrees of freedom (the real and imaginary parts). The two terms in Eq. |6.43| can 
be thought of as the separate contributions from particles and antiparticles. We 
note that \e/j,\ < ujq = V^o + m 2 here to avoid the problem of negative occupation 
numbers. For flat space with <7o = this requires \efi\ < m. 

The high temperature expansion of Eq. |6.43| presents some difficulties. There are 



a number of methods in the literature [19], |20], |2l], ^] which are quite involved. Here 
we will show in a very simple way how the Mellin-Barnes integral representation in 
Eq. |5.9| can be used to obtain the leading terms in the expansion. This method also 
has the advantage of being directly generalized to curved space as we will show in the 
next section. The first step, as for // = 0, is to expand the logarithms using Eq. p73| . 
This gives 



1 00 1 



(3^k 

^ k=l n 



V e -^n^M . ( 6 .4 6 ) 

ft k=l n ^ 



We cannot just use the Mellin-Barnes representation for the exponential in Eq. |6.46 



because the sum over k will diverge due to the presence of cosh(fc/?e/i). This means 
that we must deviate slightly from the procedure used in the non-relativistic case. 
Instead we will argue that because we are interested in (3m « 1 we also have 
(3\en\ « 1 since | e/x | < m. We therefore expand cosh(/c/3e/i) in powers of \i and keep 
only the first few terms : 

^^4EE e "^{r + ^ 2 ib + ^eVifc s } • (6-47) 

P k=l n ^ ^ 

The first term on the RHS of Eq. |6.47| is just the result for the thermodynamic 
potential found earlier for \x = (apart from an overall factor of 2). An equivalent 
way of arriving at eq. |6.47| is to start from Eq. |6.43| and expand the thermodynamic 



potential directly in powers of the small quantity (3e/i. Since fiiyo 1S an even function 
of n, only even powers of \i will arise. 

The only difference between the various terms in Eq. |6.47| is the power of k which 
enters the sum. We can define 

oo 

/i = 53E e "* /3w " A:21 " 1 ' ( 6 - 48 ) 

k=l n 
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where / is an integer. We evaluated f in Sec. |572|. Exactly the same procedure used 
there can be applied to calculate fi for / = 1, 2. Use of Eq. [5i] results in the contour 
integral representation 

fi = —. / da T(a)f3- a ( R (a + 1 - 21) ( £ , (6.49) 

^1 Jc-ioo n 

after performing the sum on k in terms of the Riemann (^-function. We have already 
evaluated Yln^n 01 i n the large box limit in Eq. |5.11| . After using this expression in 
Eq. |6.50| along with the duplication formula Eq. |5.13| for the T-function we find 



This is valid for c > max (D/2,2l— 1). 

It now only remains to study the pole structure of the integrand and evaluate 
the residues. iY 2 ^)^^ + 1 — 2/) has only a single simple pole at a = 21 coming 
from the Riemann (^-function. The residue at this pole is T{1 + 1/2). (Note that the 
potential poles of r(^2~) coimn g at a = —1, —3, —5, . . . are cancelled by the zeros 
of C,r{°- + 1 — 2Z) at these values of a. We are assuming that / > 1 here since I = 
was considered in Sec. |5T2"| .) As before the poles of r(^^) occur at a = D — 2n for 
n = 0,1,2,.... For odd D none of these poles coincide with the pole at a = 21 and 
therefore the integrand of Eq. |6.50| has only a series of simple poles. For even D we 
can have a double pole at a = 21 corresponding to n = D/2 — I when I < D/2. This 
presents enough information to evaluate // for all D. Rather than write down the 
general result, we will simply quote the results for D = 1, 2, 3. For D = 1 we have 

%^ - -\r 2 + (5- X m-\m^ - ieV) + ■ ■ ■ • (6.51) 
For D = 2 we have 

^Tyo Cr(3) 1 2 , o 2 2\ 

p — —p (m + 2e m ) 



Vr, 7T 47T 

1 

"2^ 



+ -^/T 1 (m 2 + e 2 // 2 ) ln(/?m) + • • • . (6.52) 



For D = 3 the result is 

- -^/?" 4 + ^/T 2 (m 2 - 2e 2 m 2 ) - A-/T ^(m 2 - f eV ) + • • • . (6.53) 
Vy, 45 12 071 2 

For D > 4 the leading terms in the expansion are 

5S£ * -2 7 r-^)/ 2 r(^±l)^P + l)r D - 1 (6.54) 



+ t n -(D+m T (^_±y R{D _ i)^ [ma _ (D _ l )e V] + ■ 

We will derive a more general version of this result in the next section. 
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7 Curved Space 



We will now show how to obtain the high temperature expansion for the thermody- 
namic potential of the non-relativistic and relativistic ideal gas for a general space £ 
with a possible boundary. 

7.1 Uncharged Non-relativistic Gas 
For the non-relativistic gas we have 



with u n = a n /{2m). Using Eq. |573| to expand the logarithm, and making use of the 
Mellin-Barnes integral representation Eq. |5.9| leads immediately to 



n ^ = -p^- J rf«r(a)r a C«(« + l)(^^ Q ) • (7-2) 
The remaining sum in Eq. |7.2| is just the energy ^-function of Sec 3, which was 



0J[S 



if we ignore the renormalization length. We can therefore write 

1 1 /*c+ioo 

n T ^ = --— dai>)/r Q G?(« + lV(a + l) , (7.4) 

P 2m Jc-ioo 

where we take c > D so that the integration contour is to the right of the poles of 
u(a + 1). We can analyze the pole structure of u(a + 1) exactly as in Sec. 3.3. Note 
that 



lu (a 



n n 

i r°° 

= (2m) a V— - - / dtt a - l e- ta " 
n r(a) J 

Jo 



r(a) 

where 0(£) = X]n e ~* <J ™- We therefore have 

POO 

r(a)w(a + l) = (2m) a / dtt^Qit) (7.5) 

= (2m) a /i(a) + (2m) a / 2 (a) (7.6) 
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where 



h(a) = f°dtt a - l Q(t) , (7.7) 
Jo 

poo 

I 2 (a) = / dtt^Qit) . (7.8) 

J to 

As before, to is an arbitrary positive constant. Repeating the argument given in 
Sec. |3.3| shows that l2{cx) should be an analytic function of a for any £ > 0. Therefore 
any poles of T{a)uj{a + 1) can come only from I\{a). As we did earlier, we can take 
to to be arbitrarily small and use the asymptotic expansion Eq. |3.39| for 0(t) in the 
integrand of I\{a). This gives 

oo ,a+k-D/2 

a + k — V/ 2 



=0,il, 



showing poles at a = D/2 — A; for A; = 0, |, 1, 

The only other source of poles in the integrand of Eq. |T7^ is from the simple pole 
at a = of (n(a + 1). This may coincide with the pole of T(a)uj(a + 1) at k — D/2 
if D is even, giving a double pole for the integrand of Eq. |7.4j . In any case, even if 
D is odd, we require the finite part of T(a)uj(a + 1) at a = which is impossible to 
write down in the general case. In order to avoid this calculational problem we will 
keep only the terms up to k = (D — l)/2 in the high temperature expansion of flr^o- 
We therefore find 

^-H^r D £' i£i ui+D ' 2 - k)h - <7 - io) 

fe=0,|,l,- 



Given that #o = Vs, the first term in this result is seen to agree with that of Eq. [574 
found for flat E. The remaining terms in the sum over k in Eq. |7.10| represent the 
corrections present if £ has a boundary or is curved. In particular cases they can be 
written down from the known coefficients 9k- 



7.2 Uncharged Relativistic Gas 



For the relativistic gas we still have Eq. 7A holding, but with u n = \fu n + rn 2 this 



time. The results of Sec. 3.3 are directly relevant here. If we split up u(a + 1) as in 
Eq. |3.36| and use Eq. |3.40| we find 

00 t k+(a-D)/2 

u{a) = J£ k k ^ {a _ D)/2 + analytic terms . (7.11) 
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As in the non-relativistic case, if we are to avoid terms which require a knowledge of 
the finite part of u(a + 1), we want to only keep those poles with a > 0; namely, 
k — 0, h, . . . , ^yi. This gives 

(D-l)/2 2k 

n T ^c-n-^^(3- D -i r (^- fe )u) Ul+fl-2*)0*.(7.12) 

fe=o,ii,- V J 

Again if we use 6q = Vs the first term gives a result in agreement with our flat space 
calculation. 

7.3 Charged Non-relativistic Gas 



We can begin with eq. |6.17| and expand the logarithm as usual to obtain 

00 D k0efi 



n 



^ fc=l n 

The only difference with the expression in Sec. 7.1 is the presence of e kl3etl . It is easy 
to see that in place of Eq. |0| we will get 

1 1 rc+ioo 

Q T¥0 = ~— dar(a)(3- a Li 1+a (e^)u(a + l), (7.14) 
where Li is the polylogarithm function defined in Eq. |6.21 . 



The advantage of this result over that in Eq. 7.4 is that Li\ +a is an analytic 



function of a provided that /x < 0. Therefore so long as \x < we have 

1 / \ D/2 oo / r \ k 
fe=0,i,l,... x x 

The first term of this result is in agreement with that in Eq. |6.22 . 

7.4 Charged Relativistic Gas 

For the relativistic gas we had from Eq. |6.47| and |6.48| 



fliyo * -| {/o + ^ W/x + ^/? 4 eV 4 / 2 J> • (7.16) 



The result in Eq. |6.49| still holds, so we have 



1 



c+ioo 



fi = — daT(a)(3- a ( R (a + l-2l)uj(a + l) . (7.17) 

ZTTl J c—ioo 
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The first term on the RHS of Eq. |7.1(j| is just the fi = result, and therefore should 



be twice the result found for the real field in Eq. |7.12j . The leading terms in fx and /2 



may be evaluated by using Eq. |7. 1 1| in a straightforward manner. In order to avoid 
a knowledge of the finite part of u(a+ 1) we must only keep those poles of u(at + 1) 
which lie to the right of those of (n(a + 1 — 2/); namely, those with a > 21. For fi this 
means that we can only keep pole terms with a = D — 2k for k — 0, |, 1, . . . , ^y^- 
It is therefore not possible to obtain general results if D < 3. For / 2 we must have 
D > 5. For specific spaces these restrictions can be removed; however in the general 
case we are unable to write down the general term. 
For fi assuming D > 3 we have 

fx - ^ D+1]/2 {3~ D £ | r(^±l-A;)U J D-l-2^ (7.18) 

^ TX -^m p -B T (^+ly R{D _ m + ... (7 . 19) 

if we keep only the leading order term and use #o — ^s- This last result is general 
and holds for any space E irrespective of whether there is a boundary or not, and 
independently of any boundary conditions imposed on the field. We have therefore 
established that the leading terms in the high temperature expansion of the thermo- 
dynamic potential for D > 3 are 

n*yo * -2^ D +^T(^)UD + l){3- D -^ 

lHi»W r (E^±) UD _ l)p-°[ Bl + (D- l)e 2 /.V s ] + ■ • • . 



For flat space with no boundary Q\ = —m 2 Vz and we recover the result of Eq. |6.54 



The result in Eq. ^20| was obtained in Ref. |2l| using a different method to that used 
here. 



8 Charged Interacting Gas 

We will now consider the inclusion of a self-interaction term to the ideal relativistic 
gas. The classical action will be written as 




where Kreed^l 2 ) is some potential term which represents the interaction. We will 
assume that the potential depends only on |$| 2 to ensure invariance under the U{1) 
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gauge transformation Eq. |6.24| . Following through the same steps as those which led 
to Eq. |6.35| we find 



C = $ t $ + ie/i($ t $-$ t $) + eV$ t $- |D$| 2 

- TO 2 $t$ _ c/( x )$t$ _ y trcc (|$| 2 ) , (8.2) 

as the Lagrangian density describing the theory with a conserved charge. Variation 
with respect to $ and $t results in 

= -i> + 2ie//$ + D 2 $ - (m 2 - e 2 /i 2 )$ - C/(x)$ - K rce $ , (8.3) 
= -$t _ 2 ^$t + D 2 $ f - (m 2 - eV)$ f - f/(x)$ f - V^fct , (8.4) 

where V? ree = ^ trcc (|$| 2 ). 

In order to evaluate the thermodynamic potential we will calculate the excitation 
energies [23] by looking at the split of $ into the background field $ and the quantum 



fluctuation : 

$ = $ + V- (8.5) 



This split may be substituted into the field equations |8]3| and and the result lin- 
earized in if). (Compare with the zero temperature result in Sec. 3.3.) The fluctuations 
ip obey 

= + 2ze/«/> + D 2 ^ - (™ 2 - eV + U(x))ip 

-vL e m 2 w - v:um 2 ) d*iv + *V] , (8.6) 

and the complex conjugate equation. If we assume that $ is constant as we did for 
the zero temperature theory (or less strongly that the quantum corrections to the 
derivative part of the potential are negligible in comparison with the tree level part) 
then we may proceed in much the same way as we did in Sec. 3.3. 
Expand 

V> = ]TA n e-^7n(x), (8.7) 

n 

= Y, B nZ- iEnt U*), (8-8) 

n 

where / n (x) obeys Eq. |6.38| . Here A n and B n are arbitrary expansion coefficients. 
Substitution into Eq. [81] and its complex conjugate leads to two coupled equations 
for the expansion coefficients if we use the linear independence of the / n (x) : 

= [(E n + e^ 2 -a n -m 2 ~V; vcc -m 2 V:L]A n ~^ 2 V: cc B n , (8.9) 
= [(^-e/i) 2 -a n -m 2 -K r ce-|$lXU5n-($W; ee A„. (8.10) 
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In order that these equations have a non-trivial solution for A n and B n the 2x2 



matrix of coefficients of A n and B n found from Eqs. |S.9| and |8.10| must have a vanishing 
determinant. This leads to E n = where 



Et = {Ml + eV ± [4eVM^ + |$| 4 (V t L) 2 ] 1/2 } 



1/2 



with 



Mi 



m 2 + V'+m 2 V t 



tree ' 



(8.12) 



This agrees with the results of Refs. [23|, |24] if we specialize to flat space and write 
the results in terms of real fields. As a check on the results if we take Vt ree ~~ > 0, the 
result in Eq. |3.11| reduces to those found in Eqs. |6.40| and |6.41| for the non-interacting 
theory. 

The thermodynamic potential is therefore 



n 







-PE+ 



+ In 1 - e 



(8.13) 



with E± given by Eq. |8. 1 1| . An exact evaluation of this expression is impossible for 
general S. Considerable progress can be made in flat space, but we will keep E general 
and show how to obtain the leading order behaviour at high temperature. 
Define 



Mi 



c< + SM 2 



(8.14) 



where 



6M 2 = v: rec m 2 ) + m 2 v:u\^ 



(8.15) 



and u n = \Ja n + m? as usual. By expanding E„ to first order in the potential Vt ree 
it is straightforward to show from Eq. p. 11 that 



E^ = uj n ± e/i + 



5M 2 



+ ■ 



(8.16) 



This result may be used in Eq. |8.13| with the result expanded to first order in 5M 2 . 
We find 



ttr^o = \ Yl ( ln i 1 ~ e _/3(£J "~ eAl) ] + In [l 
P 



-P(u)„+efi) 



}} 



+ Yl { [e P(tJn ~ efi) - !] _1 + [e /3( "" +eM) - 1] _1 } + • • • . (8.17) 
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The first term on the RHS of Eq. |B.17| is just the contribution from the free theory 
which was calculated in the previous section. The term in SM 2 can be evaluated in 
quite a simple manner if it is noted that 

— \e^^ - I]" 1 = 4tA> I* [1 - e-^ ±e ^l . (8.18) 
2u n 1 J 13 dm 2 1 J v ; 



If we call 



= \ ( ln I 1 ~ e' 13 ^-^] + In [1 - e -/3K+^)] } ? ( 8 . 19 ) 

n 



then 

Slr^o = + 5M 2 ^n^ + ■■■ . (8.20) 

^2yo was gi ven i n Eq. [7.20| for D > 3. If we use the fact that m 2 first enters 6i, and 
that 9^26*1 = -Ve, we have 

= V^)/ 2 r(^)c,p - + • • ■ , (8.21) 

as the leading term at high temperature. When the leading correction due to the in- 
teraction is included in the thermodynamic potential, the high temperature expansion 
becomes 

^0 * -2n-^ 2 r(^±±)c R (D + l)f3- D -iVx 
-n-^/ 2 T(D/2)( R (D)(3- D 9 1/2 

x {6 1 + V E [-5M 2 + (D - l)e 2 /i 2 ] } + •■•• (8.22) 

In the case of greatest physical interest, D — 3, if we take £/(x) = £i? with i? 
treated as constant and ignore boundary terms we have 9\ = —[m 2 + (£ — 
Dropping #1/2 (since it is a boundary term) we find 



^0 „ _^ r 4 + l r 2 



™ 2 + (£-^ + K ree (l$l 2 ) 



V33 45' 12' 

+ |<l|V t / ; cc (|<l| 2 )-2eV +•■• , (8.23) 



as the leading terms at high temperature. 
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9 Bose-Einstein Condensation for the Non-relat- 
ivistic Ideal Gas 



The fact that below a certain critical temperature an ideal gas of bosons can con- 



dense into the ground state was first predicted by Einstein |[25| . The basic physics 
of the situation can be understood using a simple argument based on the Heisenberg 
uncertainty principle. Suppose that we have N bosons confined in a volume V. Each 
boson has a volume N/V = L 3 available to it. Because the bosons are confined there 
is a zero-point energy of E = p 2 /(2m) present. The uncertainty principle relates the 
momentum associated with the zero-point motion p to the linear size of the confining 
region L by pL « h. We therefore find E ss h 2 /(2mL 2 ) = rf-(N/V) 2 / 3 . The spacing 
between the ground and the first excited state would be expected to be of the order 
E Q . If kT > E then there is sufficient thermal energy for there to be many transitions 
between the ground and the excited states; however if kT < E then we would not 
expect this to be the case. For T < T c where kT c m Eq there is insufficient thermal 
energy to excite particles from the ground state and the bosons will condense. We 
estimate that the critical temperature T c for this to happen is T c rs -^(N/V) 2 ^ 3 . 

The argument just presented is of course no substitute for a proper mathematical 
derivation; however we will show that the result obtained by this heuristic argument 
is of the correct order of magnitude. A derivation of Bose-Einstein condensation 
(BEC) along the lines followed in these lectures will be given. This has the feature 
that the connection between BEC and symmetry breaking in quantum field theory 
becomes very obvious. (The connection was described in relativistic field theory in 



Refs. pL 126 1 



The non-relativistic gas has an effective potential determined in Eq. |6.1| to be 

ficfr = Jda x {^HD^I 2 + f/(x)|^| 2 - e/z|tf | 2 j + Q T ^ , (9.1) 

with fi^yo given by eq. |6.17| . An important point is that for the ideal gas (defined 
as having no self-interactions) Qt^o does not have any dependence on ty. When we 
form the background field equation there is no contribution from Qt^o and we find 
simply 



J-D> + C /(x)-e, 



tf(x) = . (9.2) 



(We have left A and U(x) non-zero since they do not affect the main analysis given 
here.) Although we cannot solve Eq. |9.2j exactly, we can expand ^(x) in terms of the 
modes / n (x) defined in Eq. |6.14| as 

*(x) = £C n / n (x). (9.3) 
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Here C n are some expansion coefficients to be determined. Using Eq. |6.14 it is easy 

(9.4) 



to see from substituting Eq. |973| into |T2] that 

y^X^n ~ e/i)C n / n (x) = . 



Because the {/„(x)} form a complete set of functions, we must have 

= (cr„ - e\i)C n , 



(9.5) 



holding for all n. 

To avoid negative occupation numbers we had to restrict e/i < a where <t was 
the smallest eigenvalue for Eq. |6.14| . We will define a critical value of /i by /i c where 



e/i c 



a . 



(9.6) 



If /i < jj c then the only solution to Eq. |9.5| is if C n = for all values of n. This implies 
that ^ = in Eq. |9.3| vanishes and there is no symmetry breaking. 

Suppose that it is possible for /i to reach the critical value of /i c defined in Eq. 
In this case Cq in Eq. |9.5| is not determined; however, C n = for all n ^ since 
cr„ > cr in this case. The background field is given by 



*(x) = C /o(x) , 



(9.7) 



and is therefore determined by the eigenfunction of lowest eigenvalue; that is, by the 
ground state. This corresponds to symmetry breaking. 

We can make a connection with BEC by looking at the total charge, or particle 
number. From Eq. |4. 13| we have 



Q = Qo + Qi 



where 



Qo 



and 



Qi 



e /der v l^l 2 = e\C n 



on 



(9.8) 



(9.9) 



(9.10) 



This has separated off the "classical" from the "quantum" contribution. Because Cq 
is associated with a non-zero value for ^, it is natural to try to link Q Q with the 
ground state charge. 

We will now specialize to S = 13L D and set e = 1 so that we talk about the particle 
number TV rather than the charge. Again we will impose periodic boundary conditions 
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inside a box with the infinite box limit taken. a n was given in Eq. [xT] and we see 
that cr — resulting in the critical value of /x defined in Eq. being /x c — 0. The 
eigenfunction corresponding to the lowest eigenvalue, normalized to the box volume, 
is just 



/o(x) = V T 



-1/2 



(9.11) 



which is seen to be constant in this case. Thus if symmetry breaking occurs it will 
result in a constant value for the background field. 

We have already determined fi^o i n Eq. |6.20| or |6.22| . Forming iV\ = —^VLt^q 
we find 



m 



2nf3 



k=l 



k D/2 



D/2 



(9.12) 
(9.13) 



Suppose that we initially set Co = 0. Then we must have the total particle number 
given by Ni, since Nq will vanish : 



N = Vy, 



m 



2nf3 



D/2 



(9.14) 



This equation determines \x in terms of the total particle number. If // < then 
the sum in Eq. |9.12| , which defines the polylogarithm, converges for all values of D. 
However as /x — > the sum only converges for D > 2. We will assume that D > 2 
initially and examine what happens if D < 2 later. For D > 2 the sum in Eq. |9.12 
is bounded above by the value obtained when /x = 0, so that N% ~ T D / 2 for large T. 
From Eq. |9.14| this means that for high enough temperatures, no matter how large iV 
is, we can always solve Eq. |9.14j for /x with /x < 0. Thus at high enough temperature 
there is no symmetry breaking. 

Now consider what happens when the temperature is reduced from very large 
values. Because N in Eq. |9.14j is constant, as T is reduced, Li B / 2 (e /3At ) must increase. 
This means that /x must increase towards the critical value of /x = 0. Eventually T 
will reach a critical temperature T c at which /x = 0. From Eq. lO] we find 



N = Vj 



( 



mkT r 



V 2tt 



D/2 



UD/2) 



(9.15) 



or solving for T c 



2tt 
m 



xV 



Vx(r(D/2) 



2/D 



(9.16) 
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If we take D = 3 this gives the dependence on the density found from our simple 
argument at the beginning of this section. 

As T is reduced below T c , ii cannot continue to increase beyond the critical value 
of fi = 0. Thus when T < T c , ll remains frozen at the value fi = 0. It is therefore 
no longer possible to have Co = because we can no longer solve N = N\. From 
Eqs. p. 15 and |9.13| with ll = we have 

D/2 



N X = N 



Because N = N + N% we find 

N = N 

Using Eq. O, A^ = |C | 2 , we find 



rp \ D/2 



(9.17) 



(9.18) 



Co = N 1 ' 2 



D/2 



1/2 



(9.19) 



resulting in 



N 



1/2 



1 - 



rp \ D/2 



1/2 



(9.20) 



We can interpret iV as the number of particles in the ground state, and Ni as the 
number of particles in excited states. 

For D < 2 the situation is different from the one we have just described. If we 
try to let /i — > in Eq. |9.12j we find that the sum is not bounded. This means that 



regardless of how large N is, or how large T is, we can always solve Eq. |9.14| for fx 
with a value ii < 0. Physically, this means that we can accommodate any number of 
particles in excited states, and that BEC in the sense of symmetry breaking will not 
occur. It is easy to see this explicitly for D = 2 because we can perform the sum in 
Eq. |9.12| to obtain 



iVi 



-Vy, 



m 



2n[3 



ln(l 



Setting N = Ni and solving for \i gives 



kT In 



1-e 



(9.21) 



(9.22) 



Thus for any finite particle density and any non-zero temperature we can never have 
fi = 0. We always have fi < 0. For very small T we may expand the logarithm in 
Eq. to obtain 



-kTe v ^ mkT 



(9.23) 
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As T becomes smaller, /i can become arbitrarily close to /z = but never actually 
reaches this critical value. For D = 1, it is not possible to evaluate Lzi/2 in closed 
form, but divergence of the sum shows that BEC is not possible in this case. 



BEC for the D-dimensional gas has been studied in Refs. p7L 28 



10 BEC for the Relativistic Ideal Gas 

We can apply much the same formalism as we did in Sec. |9|. The effective thermody- 
namic potential is 

n eS = n<°> + n T=0 + fi T ^ , (io.i) 

where 

fi(°) = [da x {|D<5| 2 + m 2 |<5| 2 -e 2 /i 2 |$| 2 + ^(x)|i>| 2 } , (10.2) 



is the classical contribution which follows from Eq. |6.35| . f2r=o is given by Eq. |6.42 



and Qt^o is given by Eq. |6.43| . Again there is no $ dependence in Qt=o or in Qt^o 



so that the background field equation follows from Eq. |10.2| as 

= -D 2 $ + (m 2 - eV)$ + C/(x)$ . (10.3) 
If / n (x) obeys Eq. |6.38| we may expand the background field $ as 

$(x) = ^C n / ft (x). (10.4) 

n 

The result may be substituted into Eq. |10.3| and we find that 

K + m 2 - e 2 ii 2 )C„ = . (10.5) 

We will define the critical value of [i by 

e/i c = V"o + m2 • (10.6) 

For /i < /i c it is easy to see that C ra = for all n. This results in $ = and there is 
no symmetry breaking. If fi can reach the critical value /i c , then Co is undetermined 
and we have 

$(x) = C /o(x) . (10.7) 
The total charge can be decomposed as before in Eq. |9.8| where now 

Qo = 2e 2 njda x \$\ 2 = 2e 2 n\C \ 2 , (10.8) 

Qi = -§^t^- (10.9) 
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We can now specialize to flat space and take the high temperature limit, fi^yo 
was given in Sec. 6.2. For D > 3 we may use Eq. |6.54| which shows that 



\ 2 



10.10) 



is the leading term in the high temperature expansion. (In fact Eq. [7.20| shows that 
this is true even for a curved space, provided that we are justified in dropping the 
zero temperature term.) From Eq. |10.10| we see that for large T, Qi « T D ~ l . This 
means that for high enough T we can always solve Q — Q\ for fi with \i < \i c . In 
this case $ = and there is no symmetry breaking. As T decreases \i must increase 
if Q is to remain fixed. However /i can never increase beyond /i = /i c . Thus we may 
define a critical temperature T c by 



q ~ ^ D +^V(^±±)UD - l)(fcT c )^V s 2eV 



or explicitly, 



Q 



e 2 /j, c Vza D 



1/(D-1) 



where 



For T < T r we have 



np =2K HD + Q/* r (E+l) UD _ 1)m 



D-l 



10.11) 



10.12) 



(10.13) 



(10.14) 



from Eq. |10.10| since fi = fi c for this temperature range. Because Q = Q + Qi we 
have 



Qo = Q 



D-l 



10.15) 



Using Eq. |10.15| in Eq. |10.8| with /j, — fj, c shows that 

1/2 



C = (\<*dV^ [(kT c ) D ~ l - (kT) D - 1 ] . 



(10.16) 



For flat space a = and e/i c = m from Eq. |10.6| . Also /o( x ) = > so that 

1/2 



10.17) 
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gives the background field and symmetry breaking can occur. 

We can analyze the flat space situation in another way which does not involve the 
high temperature limit, and which provides a suggestion for the general approach we 
will discuss later. Beginning with 

n T ^ = - ( ln i 1 ~ e"^-^] + In [1 - e -/»( w »+ e <*)] } , (10.18) 

n 

and expanding the logarithms gives 

1 00 1 

Q T ^ Q = -±.J2j; (e fc/3e/i + e~ k ^) ^ e~ k ^ . (10.19) 

" k=l n 

In the large box limit we established in Eq. |5.11| that 

J2 = V^j e -*th0W (1Q 20) 

n J \ ) 

Although the integral cannot be done in terms of elementary functions, it can be done 
in terms of Bessel functions with the result 

/ \ (£>+l)/2 

jry*** = 2(3kV s (^-j K {D+1)/2 (kf3m) . (10.21) 

We therefore obtain 

/ \ (D+l)/2 oo 

flzyo = ~2Vx g AT^ +1 >/ 2 (e k ^ + e~ k ^) K (D+1)/2 (k(3m) .(10.22) 



No approximations have been made to obtain this expression other than the large 
box limit. Although this expression is not very well suited to examine the high 
temperature limit (although it can be done f29|), it can be used to see the restriction 
on the allowed dimensions for BEC. 
Using Qi = -^fijyo we find 

/ \ (£>+l)/2 oo 

Ql = 2ePV * ^ k ~ iD ~ 1)/2 - e '^) K {D+m (kf3m) . (10.23) 

Provided that Q\ remains finite as /i — > /i c then our analysis of BEC presented above 
is alright. Convergence of Qi is reduced to a study of the sum over k, which is 
determined by the large k behaviour of the summand. For large enough k we can use 
the asymptotic expansion for the Bessel function, 

K„(z) ~ ^ e- z {1 + Oiz- 1 )} . (10.24) 
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The convergence of the sum is determined by the behaviour 

/ \ (D+l)/2 oo , 

«■<" - - fe) S "-^^^ ■ (io - 25) 

where terms which are finite as fi — > /i c or which are less divergent than the one 
indicated have been dropped. Thus, 

- A*c) - eV s (^) £ k-DWtom) , (10.26) 

determines convergence. It is now easy to see that this sum only converges as e/x — > 
e/i c = m for D > 3. For £) = 1, 2 it does not converge. This means that for D — 1, 2 
it is possible to put any amount of the total charge in the excited states and BEC in 
the sense of symmetry breaking will not occur. The discussion, and in particular the 
restriction on the spatial dimension, is the same as in the non-relativistic case. 

11 BEC for Nonrelativistic Bosons in a Constant 
Magnetic Field 

We will now study the behaviour of relativistic and non-relativistic charged spin-0 
fields in a constant magnetic field in a flat space, £ = M D , of general dimension. We 
are interested in how the presence of an externally applied field affects the results 
found for free fields. 

The first point worth making concerns what is meant by a constant magnetic 
field in a space of general dimension. The magnetic field in D spatial dimensions is 
described by a constant antisymmetric tensor F^, not in general by a vector. (The 
electric field components F 0i are described by a vector.) In the special case D = 3, 
given Fij we can define a vector with components B l by F^ = €ijkB k where e^fc is 
the Levi-Civita tensor, (eyjt = +1(— 1) for ijk an even (odd) cyclic permutation of 
123; e^k = otherwise.) If D = 3 we are always free to choose our coordinate axes so 
that the vector B lies along one of the axes. There is therefore no loss of generality 
by assuming only F 12 = —F 2 i to be non-zero. For general D the situation is not so 
simple. 

One easy way to see that for general D > 3 it is never possible to find coordinates 
such that only F 12 = —F 2 x is nonzero is to consider the invariants which can be formed 
from Fij. For example when D = 4 there are two independent invariants :F^F^ 
and eijkiF 13 F kl . (See Ref. [[3(J for example.) If we keep only one component of the 



magnetic field non-zero then the second invariant vanishes. It should be clear that for 
D = 4 we must keep two independent components of Fy non-zero, say F\ 2 = —F 2 \ 
and F u = -F 43 . 
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For general D, rather than consider the invariants directly, it is easier to consider 
reducing F^ to canonical block diagonal form. Consider first the case of D even. 
Define D = 25 for some integer 5 = 1,2,.... We can always put F^ into the form 

Fij = diag(Ai J, . . . , X s J) , (11.1) 

where 

J = ( -i o ) ' (1L2) 

and ±zAj, . . . , ±i\s are the eigenvalues of F^. There are therefore 5 independent 
invariants we can construct from F^. We lose no generality in setting all but 5 
independent components F± 2 , F34, . . . , F2S-128 to zero. Define 

Bj = F 2j - 12 j , (11.3) 

for j = 1,2,... ,5. Then £>i, . . . ,B$ represent the independent components of the 
magnetic field. We can choose a gauge such that the non-zero components of the 
vector potential are 

4y-i = -B jX 2j , (11.4) 

for j = 1, . . . , 5. 

For D odd we may define D = 25 + 1 for 5=1,2,.... (D = 1 does not have any 
magnetic components for Fij.) In this case det Fy = since F^ is antisymmetric. 
Thus there is always at least one zero eigenvalue. The subspace spanned by eigen- 
vectors orthogonal to the one with zero eigenvalue is 25-dimensional. Therefore as in 
the even dimensional case there are 5 independent invariants. We may still choose 
the gauge where Eq. |11.4j holds. 



In order to study how the magnetic field affects BEC, we will keep only the first p 
components in Eq. |11.4| nonzero. We require p < 5. With the gauge choice in Eq. |11.4 



we require the eigenf unctions and eigenvalues of 

" - ~L t [foH. + + 4] -±±0l (11.5) 

j=l j=2p+l 

D = 25 or 25 + 1 depending upon whether D is even or odd. The advantage of our 
gauge choice is that the first p terms on the RHS of Eq. |11.5| are equivalent to a 



sum of p non-interacting simple harmonic oscillator Hamiltonians. If we consider a 
single term, say j, in the sum, then the eigenvalue is (2nj + l)eBj where rij = 0, 1, . . . 
with degeneracy eBjL 2 j_iL 2 j / {2n) . (See Ref. |yj for example.) Here we are again 
imposing box normalization with Lj the length of the box in the j th direction. If we 
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impose periodic boundary conditions in the remaining directions, then the eigenvalues 



of — 7T-D 2 are 

2m 



j=l j=2p+l V 



;n.6) 



where fc,- = 0, ±1, ... , with degeneracy 



n e -^ L ^-^ 

3=1 



11.7) 



The thermodynamic potential involves 



oo oo 







D-2p 



In [1 



-/3(<*>nk-e/j)l 



11.81 



if we take the large box limit. In Eq. |11.8| we have 

Wnk = ^(2n, + 1)^ + -. 



3=1 



11.9) 



Expanding the logarithm in Eq. |11.8| as usual, and performing the integration over k 
results in 



n 



Vy / Aefi 



oo oo oo 



-zyo 



i n? E -EE 



\- LJ - 2tt , 

\J=1 / ni=0 n p - 



' ^ J (2tt) 

Z=l ^ v ; 



)D-2p 



Vy I T-r eB; \ f m 



ft 11 9tt 



(D-2p)/2 oo 



8 \ A1 2vr / V2tt/3 

\3=1 / 

oo oo 

E--E e ~^ E '- i<2 " j+1,eB 



,lf3ep, 



X 



(11.10) 



ni=0 n„=0 



It is easy to perform the sums over m, . . . ,n p since they just involve geometric series 
and obtain 



Vs f m \ 



(D-2p)/2 oo 



where we have used 



■D/2J/3e(p-^ c ) 



2tt 



i=i 



e/i c = o"oo 



3=1 



1-e" 



11.11) 



^ 2m 



3=1 
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Defining Q 1 = — J^fijyo as before, we have 



Qi = eV s 



m 



(£>-2p)/2 oo 



27T(3 J 



D/2JPe(p-iJt c ) 



2tt 



-1 



;ii.i3) 



We can now study whether or not BEC can occur. 

From our earlier discussion, in order to have BEC we must have Qi bounded as 
ji — > /x c . Since /3eBj/m > is assumed, we may use the inequality 



1 < (l-e-T 1 < fl-e-*)- 1 



for 1 = 1,2, 



to see that 



Qi > eV s 



_m_\ (D ~ 2p)/2 
2^J 



nH)E"- 



D/2 e J/3e( M -Atc) 



11.14) 



11.15) 



This shows that Qi is not bounded as /x — > /x c for y — p < 1. When this inequality 
holds BEC does not occur because any amount of the total charge can be accommo- 
dated by setting Q = Qi and solving for /x with /x < /x c . For D = 25 this will be 
the case if p = 8 or 8 — 1. For D = 25 + 1 this will only be the case if p = 5. This 
shows that for a general magnetic field in a space of any dimension BEC in the sense 
of symmetry breaking will not occur. For special values of Fij BEC is still possible 
provided that the inequality y — p < 1 is violated. 

In order to see that BEC is still possible for | - p > 1 we need to prove that Q\ 
in Eq. |11.13| remains bounded as /x — ► /x c . This follows easily by using the inequality 
in Eq. |11.14| again to show 



Qi < eV^ 



m 



2n[3 



(D-2p)/2 




dcD, 



-1 



OO 

E 

1=1 



lP-D/2 e iPe(„-„ c ) _ ( 1L1 g) 



The sum on the RHS is bounded by ( R (D/2 - p) if f -p > 1. If D = 28, BEC can 
occur for p = 1, 2, . . . , 8 - 2 provided that 8 > 3 (ie D > 6). If D = 28 + 1 BEC 
can occur for p = 1, 2, . . . , 8 — 1 provided that 5 > 2 (ie D > 5). If we specialize to 
the case of a single component field, p = 1, our results show that in order for BEC 
to occur we require D > 5. In particular BEC in an externally applied field does 
not occur for D = 3. It is easy to show that letting Bj — >• in Eq. |11.13j reduces 
the expression to the result found earlier for the free Bose gas. There is therefore a 
substantial change in the system caused by turning on a constant magnetic field, no 
matter how small. 



Suppose now that the conditions for BEC are met 
temperature T c is defined by 



ie 



D 



p > 1). The critical 



Q = eV 11 



( mkT c 
V 2tt 



(D-2p)/2 oo 

i=i 
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n eB J 
2tt 
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11.17) 
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Here (3 C = (/cT c ) _1 . It is not possible to solve this expression explicitly for T c as we 
did for the free Bose gas. However, it is possible to obtain approximate results in the 
limits of strong and weak fields. 

Assume first that (3 c eBj » m. As a first rough approximation we can neglect 
the exponential term in Eq. [11.171 and obtain 



mkT c 
2tt 



(D-2p)/2 



D 



3=1 



2tt 



11.18) 



which is easily solved for T c . If T < T c then using the approximation in Eq. |1 1.18 
have 



we 



Qi 
Qo 



rp \ D/2-p 



Q 



D/2-p 



11.19) 
11.20) 



Comparison of these results with those in Eqs. |9.17| and |9.18| obtained for Fij = 



shows that the presence of the magnetic field has the same effect as lowering the 
spatial dimension from D to D — 2p. If we define D e ^ = D — 2p, then the condition 
for BEC is D c s > 3 for p = 1, exactly as in the free Bose gas. Physically the fact that 
a non-zero magnetic field reduces the spatial dimension by two is clear if you keep in 
mind what happens classically. For D = 3 the classical motion of a charged particle 
is a spiral around the direction specified by the magnetic field. The only free motion 
of the particle is in the direction of the magnetic field; the motion perpendicular to 



the magnetic field is confined. We will return to this in Sec. [14 



Next assume that /3 c eBj « m. A naive expansion of the exponential in Eq. |11.17 
will not work because no matter how small f3 c eBj/m is I eventually becomes large 
enough so that the sum over I diverges. In Ref. |32[] we approximated the result by 
replacing the factor of [1 — exp(— lf3 c eBj/m)]~ 1 with its upper limit obtained when 
1 = 1. This allowed an analytic result for T c to be obtained. Here we wish to show 
how a more accurate analysis can be used to obtain a better approximation. (This 
was suggested to me several years ago by Klaus Kirsten.) 

We will begin with the expression for Qi found from Qi 
used for Qt^o- We will also define 



-|l^o with Eq. pTT0| 



m 



(11.21) 



In addition, because we want to let /i — >• ji c where /i c was defined in Eq. p. 1.12 
will let 



we 



H = Hc{\ - e) . 



11.22) 
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In terms of Xj we have 



1 P 

fe= 5 ^. (11-23) 



2 



It is easy to show that 



D/2 

m x 



M/ \j"=l / ni=0 n p =0 1=1 



P 

e 

x 



exp-Z^(n J - + -)x 3 - , (11.24) 
when written in terms of the dimensionless variables Xj and e. We can now use the 



Mellin-Barnes representation in Eq. |5.9| for the exponential and perform the sum on I 
in terms of the Riemann (^-function. The sums over m, . . . ,n p can be done in terms 
of a multidimensional generalization of the Hurwitz ^-function called the Barnes £- 
function [33] . We will first outline very briefly the properties needed for our analysis 
here. 

The Barnes ^-function is defined by 

oo oo 

( B (s, a|x) = J>-X> + n. x)~ s , (11.25) 

ni=0 n p =0 

where x and n are p-dimensional vectors and a > 0. When p — 1 the Barnes (- 
function reduces to the Hurwitz ^-function. Using the identity in Eq. p,14| we can 
easily perform the sums over n and obtain the integral representation 

Cb(Mx) = — t / dtt s - l e- at X\{l-e- tx >y l . (11.26) 

All that we require are the poles of (b{s, a|x). These can be obtained in a simple way 
by following the same procedure as we used in analyzing Eq. |3.34| . The integration 
range can be split up into two parts; the first extending from to to and the second 
from to to oo. It can be argued that all of the poles come from the integral from to 
to, and as before we only require the behaviour of the integrand for small t. It is an 
easy exercise to show that the first few terms are given by 



II (1 - e"^) 



-1 




+ •■• } . (11.27) 
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It is now straightforward to show that Cb( s > Q | x ) has simple poles at s = p,p — l,p 
2, . . . . The residues at the first few poles are 



Res (b(s, a|x) 



m n 



'11.28) 



Res Cb(s, a|x) 



Yl x 3 - 2a 



= P -i 



2r(p - 1) n 
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Res Cb(s, a|x) 
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11.30) 



Clearly we must have p > 1 for the result in Eq. |1 1 . 29| and p > 2 for the result in 
Eq. [IT301 . 

From Eq. |11.24| and the definition of the Barnes (^-function, we find 



/ m \ D/2 [yj \ 1 f c+io ° 



D, 
2 ' 



3=1 



11.31) 



Closing the contour in the left hand plane and evaluating the residues using the results 
just obtained for the Barnes ^-function we find 

where we have kept only the first order correction in the magnetic field. (Recall that 
we are assuming that BEC does occur, which requires D > 5, so that the Riemann £ 
functions are finite.) We can now let e — > and set Qi = Q which defines the critical 
temperature T c : 



Q^eVx 



mkT c 
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If we let Bj — > 0, it is clear that T c reduces the result found earlier in the absence of 
a magnetic field. (See Eq. |9.15| for the particle number for the free field case.) 

In order to see more clearly how the presence of a magnetic field alters the result 
from that for the free field Bose gas, let To be the critical temperature in the absence 
of a magnetic field. (This was the temperature called T c in Sec. ||.) It is easily seen 



from Eq. 11.33 that 



T c - T - 1 C ^: 1} y(^). (11.34) 

This result shows that the presence of a weak magnetic field lowers the critical tem- 
perature with respect to that found in the absence of a magnetic field. This conclusion 



was also reached in Ref. [32] using a less refined method which led to a cruder ap- 



proximation for T c in terms of Tq. 



One final comment we wish to make concerns the nature of the Bj — > limit. The 



steps culminating in Eq. p. 1.34] assume that D is at least 5. When D < 5 we have 



seen that BEC does not occur when a magnetic field is present. This means that the 
physically interesting case of D = 3 is qualitatively different when a magnetic field 
is present (no BEC) and when a magnetic field is absent (BEC). The existence of a 
magnetic field, no matter how small it might be, destroys the BEC of an otherwise 
free gas. This was pointed out originally by Schafroth [^] who noted that it is not 
possible to treat the constant magnetic field by perturbing around the zero field result 
in any finite order of perturbation theory. It would be a very interesting problem to 
study the behaviour of a system of bosons when a magnetic field is turned on or off. 
The restriction on the spatial dimension D > 5 for the single component magnetic 
field was first noted by May MB . The general field was first presented in Ref. |32|, B9] . 



12 Meissner-Ochsenfeld Effect 

It is a well-known feature of a superconductor that when a magnetic field is applied 
to a sample there may be no field inside the sample if the field is larger than a certain 
critical value. A beautiful analysis of the Meissner-Ochsenfeld effect and its rela- 



tionship with BEC was presented by Schafroth |[34|| . Our aim here is to show how 
the approach we have been using in these lectures may be used to study this prob- 
lem. There has been some confusion in the literature following Schafroth's (correct) 
analysis. One advantage of adopting the effective action approach is that there is no 
possible source of error concerning what Schafroth called the acting field and what he 
called the microscopic field. Another source of error, which seems to have originated 



with Ref. [|35] concerns the choice of units. We will adopt Heaviside-Lorentz ratio- 
nalized units as is conventional in quantum field theory. (A discussion of the various 
units and their inter-relationship in this problem is contained in Ref. [[37].) 
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The complete thermodynamic potential is (see Sec. []) 

fieff = 0™ + ^ (0) + fir^o , (12.1) 

where 

fiem = [da x [~F i:i F ij - J^A} , (12-2) 

has been ignored until now because it is independent of ^. We will be interested in 
the case where BEC does not occur, since this is the case for the physically interesting 
situation of D = 3. The Maxwell equation, which follows from setting the variation 
of fi e ff with respect to A, to zero is 

V;i^' = <4t + 4d , (12-3) 

where 

«)r#o = -Jda x J* nd 6A t , (12.4) 

defines the current density induced by quantum corrections to the classical theory. 
(If ^ then there would also be a term from f^°) on the RHS of Eq. |12.3| .) 

We have computed Q as a function of rather than Ai, so it will be best if we 
can express J? nd in terms of a derivative with respect to Fij rather than Aj. If we vary 
Qt^o with respect to F^ we have 

Sn T ^ = [da x 5 -^- 5F l3 (12.5) 



jda x 6 -^(v i 5A j -V J 6A l 
2 [da x V j ( 6 -^]5A i . (12.6) 



Comparison of Eq. |12.6| with |12.4| shows that we can write 



J', d = -2V, ( ^ ) . (12.7) 



The Maxwell equation |12.3| becomes 

Vjfl* = 4 t , (12.8) 

with 

AO 

#y = + . (12.9) 

5Fy 
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This is a generalization of the usual vector electromagnetic result of H = B — M 
in D = 3 where M is the magnetization. In fact for D = 3 it is easy to show that 
Eq. |12.9| reduces to this standard result |37| . 



The factor of 2 which occurs in Eq. |12.9| is because we have treated all of the 
components of Fij as independent when performing the variation. However when we 
actually computed flr^o we used F 2 \ = —F 12 and defined F 12 = Bi, and so on for the 
other independent components. If we define H u = Hi and so forth in an analogous 
way, then we obtain 

H i = B i + S -^ i (12.10) 

where the factor of 2 no longer appears. (The factor of 2 can also be seen to disappear 
between Eq. |12.9| and |12.10| if we note that Qt^o only depends on through F^F % ^ 



and use FijF ij = 2(Bf + ■ ••£*).) 

We will define the magnetization by 

5Q 



Mi = ■ (12.11) 

oBi 



From Eq. |12.8| it is seen that H l i or Hi is the field set up by the external current, and 



therefore should be interpreted as the externally applied field. F %] or Bi represents 
the field present in the sample which is felt by the particles. These two fields are 
of course different because the externally applied field causes the charged particles 
to move, which gives rise to an induced current J- md in the sample. This induced 
current changes the overall magnetic field. For a Meissner-Ochsenfeld effect to occur 
we need to know if it is possible for Bi — > when Hi is changed. If this can happen 
it corresponds to the particles in the sample feeling a vanishing magnetic field. 

We will now concentrate on the case of D = 3, and evaluate the magnetization 
Mj. From Eq. |11.12| we have 



as the critical value of the chemical potential. We know from our analysis in Sec. [TTJ 
that BEC does not occur, and therefore that /1 can never reach the critical value of 
/i c . We will define 

// = /i c (l-e), (12.13) 

for dimensionless parameter e with e > 0. In order to calculate the magnetization we 
will first obtain the Mellin-Barnes type integral representation for Qt^q, similarly to 



the representation found for the charge in Eq. |11.31| . From Eq. |11.10| we have (with 
p — 1 and D = 3) 
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This may be written in terms of e, and x = (3eB/m defined in Eq. 



PT2II 



We have 
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Rather than just apply the Mellin-Barnes integral representation directly to the 
exponential, it is better to separate off the n = term for special treatment. The 
result can be expressed in terms of the polylogarithm, 
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where the first term has arisen from the n = term in Eq. |12.15| . Now we will use 
Eq. |5.9| to convert the double sum into a contour integral with the result 
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(We take c > 1 here.) The advantage of separating off the n = term as we have done 
is that we ended up with Cff(a, 1 + f ) rather than |). This has the practical 

advantage that we can let e — > without worry in all terms arising from the contour 
integral since they are well-behaved in this limit. More care must be exercised if the 
result is written in terms of Cff( a ; f )• In addition we can interpret the two terms in 
Eq. |12.171 physically. The first piece involving Li3/2 comes from the n = term and 
therefore is associated with the ground state. The remaining part comes from n > 1 
terms and is therefore the excited states contribution. 

We now close the contour in the familiar way to obtain the asymptotic expansion 
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(It is simple to write down more terms in the expansion if they are needed.) The 
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charge density is 
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(The differentiation is performed with i? held fixed.) The magnetization is 
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m dx B de 
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since // is held fixed. (It is easy to compute jrk holding \i fixed if we use Eqs. |1 2 . 1 2 



and 12.13.) It is found that 
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(12.21) 
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Now consider what happens as T is reduced from large values towards zero. If BEC 
could occur then this would correspond to e — > because \i — > // c would be reached 
at some critical temperature. However e can never reach zero because Li 1 / 2 (e~2 <E:r ) 
which enters the charge density diverges in this limit, and we saw that this is what 
prohibits BEC as symmetry breaking from occurring. Suppose that we define po to 
be the charge density in the absence of a magnetic field. Then from Eq. |9.15| we have 



Po 



in 



3/2 



27T/5 



Cb(3/2) 
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Furthermore let To be the critical temperature in the absence of a magnetic field. 
Then from Eq. |9.17| we have 



rp \ 3/2 



12.23) 



We expect that asT^ T c , e should start to get very small in the presence of a weak 
magnetic field. We can use 



Lii/2 (e 2 



'27T 

'— + ( R (l/2) + 0(ex). 
ex 



12.24) 
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This is easy to establish using the Mellin-Barnes expression in Eq. 53 as follows. 
Begin with 



Xl e~ n0 



n=l 

= — / daT (a )Cr (a+- ) 6~ 

2m y c _ ioo V 2) 



+ &(l/2) - Cb(-1/2)0 + 



We can therefore approximate Eq |12.19| by 



p ~ p + e 
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2ttx 
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Using the same expansion for the magnetization results in 

3/2 



M 



2ttx 



2m V27T/?, 

for small e. We can eliminate e between Eq. |12.26| and Eq. |12.27] to obtain 
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for T < T c . Finally we may use this result in eqs. |12.10| and |12.11| to see that 

3/2" 



H = B + ^ 
2m 



-(0 
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This shows that once T reaches the temperature T c , which is the critical temperature 
for BEC in the absence of a magnetic field, it is possible to have B = for H > H c 
where 



2m 



3/2' 
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This corresponds to an expulsion of the applied external magnetic field from the 
sample, which is the Meissner-Ochsenfeld effect. 
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13 Relativistic Charged Bosons in the Presence of 
a Constant Magnetic Field 

The starting point is the thermodynamic potential fi^yo expressed by 



13.1 



For relativistic charged bosons in a constant magnetic field it is easy to show that, 

(13.2) 



using the same notation as in Sec. |TT 

^"nk = V Onk + 



where 



a nk = ±(2n j + l)eB j+ £ (^Y . 

j=l j=2p+l ^ J ' ' 

If we expand the logarithms in Eq. |13.1| and take the large box limit we have 

M i=l \9=1 / ni=0 n„=0 



13.3) 
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|D-2p 



exp I — i/3 



k 2 + m 2 + y^(2nj + l)e-B,- 
j'=i 



l/2> 



. (13.4) 



(Compare with eq. |10.19| for the free Bose gas.) We can now use Eqs. |10.20| and |10.21 
to evaluate the integral over k : 
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We can compute Qi = — ^f2jy to be 



Qi 



V T ,e(27c) {D ~ 2p+1)/2 f3-( D - 2 P- 1 )/ 2 

x UPen _ e -lf3e^ i~(D-2p-l)/2 
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^(2^- + l)eSj 

p 1 V2 N 

+ J2(2 nj + \)eBj 
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(D-2p+l)/4 
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To examine what happens as \x — > /i c we can use the same analysis as we used at 
the end of Sec. [TD. All we need to know is the behaviour of the sums which occur in 
Eq. |13.6| . For large I we can use Eq. |10.24j and note that for all terms with rij ^ the 
sum will converge exponentially fast. This means that the behaviour of Q\ as fi — > \i c 
is determined by the n\ = ■ ■ ■ = n p = term in the sums. Thus 



(£>-2p)/4 



lP-D/^ e -We(iJ,c-n) _ (13.7) 



Z=l 



This expression only converges in the limit /i = /i c for D/2 — p > 1. This is identical 
to the result found for the nonrelativistic Bose gas in Sec. II. In particular, we can 
define an effective spatial dimension by D c q = D — 2p. The requirement that BEC 
occurs requires D c{i > 3. For p = this reduces to the usual one for the free Bose 



gas. An analysis of the Meissner-Ochsenfeld effect can be given [38, 39, 40, 41, 37 



14 



General Criterion for BEC Interpreted as Sym- 
metry Breaking 



We have seen in Sec. [TT| and Sec. [H| how the presence of a constant magnetic field 
with p independent components alters the criterion for BEC from that found for the 
free Bose gas. If we define an effective dimension D c g = D — 2p then BEC occurs only 
for D c g > 3. In this section we will see in a fairly simple way the general features 
present which alter the effective spatial dimension from D to -D e ff- This analysis was 



first presented in Refs. [37, 42 . 



14.1 Non-relativistic Gas 

Suppose that we consider any system for which the energy levels can be expressed 
as the sum of a discrete part, which we will denote by E^, and a continuous part 
which can be dealt with by box normalization. (A more general possibility for the 
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continuous part can be considered by using a method we will describe later.) We will 
take the infinite box limit at the end. Assume that the box has dimension q and write 



p just denotes the set of labels for the discrete part of the energy spectrum. The 
thermodynamic potential is 

^Vo = i^ln[l-e-^—)] 



_V q ^ f d q n ^ 1 



Kv^ f d q n I e _ kp{E d >+ g__ efl) 



3 J (2tt)« ^ k 
' P J v I k=i 

T/ / \ 3/ 2 00 

( ™ 1 V V k-^e-^^ , ( 14.21 



EE*- 

7 p fc=i 



after expanding the logarithm and performing the integration over n. The critical 
value of \i is given by 

efi c = E d , (14.3) 

where Eq is the smallest energy eigenvalue in the discrete part of the spectrum. 
Because the behaviour as // — > /x c is of interest we can separate off the lowest mode 
contribution in Eq. [L4.2| and write 

VLtjLq = fifyo + ^zyo > (14.4) 

with fify the contribution coming from p = 0, and fi^y the contribution coming 
from p ^ 0. Physically fify contains the ground state contribution and fi^yo con " 
tains the contributions from excited states. If c?o is the ground state degeneracy we 
have 

OSyo = 4. (j^) ^ U 1+q/2 (e-*0o-» >) . (14.5) 

We also have 

g/2 



1 ^ fc=l p/0 



where the sum over p is restricted to excited states only. The advantage of separating 
off the lowest mode as we have done is that even if we let // — > \i c the argument of 
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the exponential in Eq. 14.6 remains negative ensuring convergence of the sum over k. 
This can be seen by noting that for large k we have 



5> 



-k/3{Ef-en) 



where Ef represents the first excited energy level. It follows that for large k the 
summand behaves like 

k -l-q/2 e -kf3{Ef-e^) _ 

As e/i — > e\i c = Eq the sum over k is exponentially damped because E\ — Eft > by 
assumption. This exponential damping persists even for fi = /i c , and is also true even 
if we differentiate with respect to /i before taking the /i — > fi c limit to get the charge. 
Thus, whether or not a phase transition representing BEC occurs is determined by 



We can write 



With Qf = - Jjfif^o we have 



Qi = Qf + Ql 



Qf = eV q d 
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(Q\ x can be defined in terms of fi^y in an obvious manner.) We have argued that 
the behaviour of Qi as \i — > \i c is determined by Qf , so 



Qi(fjL fic) - eV q do 



m 



27T/3 



q/2 



Li q/2 le-P<^ 
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Li g /2^e pe ^ c ^ j is only finite as \i — > \i c if q > 2. For q — 0, 1,2 Q\ will diverge as 

ji — ► n c meaning that BEC cannot occur as it does for the free Bose gas. The actual 
divergence as \i — > \i c is easily obtained from the behaviour of the polylogarithm [07 



For a magnetic field with p independent components in D spatial dimensions we 
had q = D — 2p as the dimension associated with the continuous part of the spectrum. 
(See Eq. |11.6 in the large box limit.) The restriction q > 3 found here is therefore 



fully consistent with our earlier result in Sec. |TT] for the special case of a constant 
magnetic field. 



14.2 Relativistic Gas 
In this case we have 



n 



5> 



-0E+ 



+ In 
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-PE- 
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where 

= V<?n + m 2 ± e/i . (14.11) 

Assume, as in Eq. |14.1| , that a n splits into the sum of a discrete part and a 
continuous part which we deal with by box normalization : 



On = cr p 



1=1 



If we expand the logarithms in Eq. 14. 1C and take the infinite box limit we obtain 



n T ^ = -YjL ^ J2 I (e kPe » + e~ k ^) J _|^e- fe/ V CT °+ n2 + m2 . (14.13) 



We can do the integration over n in terms of a Bessel function as in Sec. 1C with the 
result 

A = _l|(2 7r /3)( 1 -'')/ 2 yA;-( 1+ ^ 2 y(^ + m 2 )fe +1 )/ 4 

r k=l p 

xK {q+1)/2 (kPy/<r* + m 2 ) (e~ k ^ + e kl3e ^ . (14.14) 
As in the nonrelativistic case we will separate off the p = term and define 

= _M( 27r/ 3)( 1 ^)/ 2 ^A;-( 1+ ^ 2 (^ + m 2 )^ +1 '/ 4 

P k=l 

xK {q+1)/2 (V^ + m 2 ) (e~ k ^ + e k ^) . (14.15) 

fi^Yo is defined as in Eq. |14.14| , but with the sum restricted to p ^ 0. Due to 
the behaviour of the Bessel function for large values of its argument, fi^ remains 
finite (as does its derivative with respect to /i) even for /i = fi c . The convergence or 
divergence of the charge Qi is determined by the behaviour of ^fy - 
For large k we have 

/ \ V 2 

K( q+ i)/2 (k(3^v d + m 2 ) ~ {a d + m 2 )-VV^V^ . (14.16) 

Noting that e/i c = a/o"q + m 2 we have the large behaviour of the thermodynamic 
potential determined by 

nf^-AO - -^(^) ?/2 f:^ 1 - 9/2 e-^-^ (14.17) 

3/2 













( e/i c 


/3 





J Li 1+ff/2 (e-^-^J . (14.18) 



Just as in the nonrelativistic case, it is easily seen that Q\ remains finite as fi — > fi c 
only for q > 2. 
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15 BEC for the Self-interacting Theory 

We will assume that D = 3 here, and that BEC does occur. (This requires a contin- 
uous energy spectrum as described in the last section.) We will also assume that $ 
is constant. All of these conditions are met for S = R 3 . For the self-interaction we 
will take 

i/ tree (|<i| 2 ) = ^|<i| 4 , (i5.i) 

6 



to be the usual \<p 4 theory. Using the results of Sec. |373| with the zero-point energy 
term dropped as irrelevant at high temperature, we have the total effective thermo- 
dynamic potential given by 



Vx v * r " 1 6' 1 45 

+i/T 2 [m 2 + (t - \)R - 2eV + y |$| 2 ] + ■ ■ ■ • (15.2) 

The background field $ must satisfy ~^^es = 0, giving 

~ (m 2 + £R - eV)^ + ^l$| 2 $ + 4/?~ 2 $ • (15.3) 

3 18 

There are two possible solutions to this : $ = 0, corresponding to the symmetric 
phase, and 

\*\ 2 = \{#f-m 2 -iR-±fr*), (15.4) 
which corresponds to the broken symmetry. The total charge density is 

Ve op 

~ 2e 2 /i|$| 2 + -e 2 p/3~ 2 (15.5) 
3 

If we define 

po = 2e 2 /i|$| 2 , (15.6) 
Pi = \e 2 p(3~ 2 , (15.7) 

we can interpret po a s the charge density of the Bose condensate. In the symmetric 
phase we have <£> = so that 

p ~ l -e 2 p(3~ 2 . (15.8) 
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As T decreases from large values, for fixed total charge /x will increase until we reach 
the critical value set by 



e [i 



(15.9) 



at the critical temperature T c . This is the temperature at which the symmetry is 
broken. For T < T c we have (see Eq. |15.4 ) 



l$l 



(15.10) 

and \i = /i c . Setting T = T c in Eq. [15. 8| we find 

e/i c = % c 2 . (15.11) 

Using this in the left hand side of Eq. |15.9| gives us an equation which completely 
determines T c in terms of the charge density. From Eq. |15.7| we have 



Pi ^ P 



(15.12) 



and hence 



Po- P 



15.13) 



This is exactly the result found for the free Bose gas (although of course the result 
for T c is different). 

It is possible to study the nonrelativistic limit of the theory. This has been done 
in a very beautiful paper |24]. It is also possible to generalize to a non-abelian gauge 
symmetry [E3 . 



16 BEC in an Isotropic Harmonic Oscillator Po- 
tential 

We will now consider the case of non-relativistic bosons in a harmonic oscillator 
potential. This simple theoretical model provides a good description of the magnetic 
traps used in the recent experiments on BEC [[0], ^5], fl6|, . Interactions among the 



atoms will be neglected. As a first approximation this would be expected to be not 
too bad because the atomic gases are dilute. The grand canonical ensemble will be 
used. 
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We will only look at D = 3, although the analysis may be generalized to other 
spatial dimensions. We will take 

£/(x) = \m(u\x 2 + uly 2 + u 2 3 z 2 ) , (16.1) 

where u)\ , uo 2 and uo 3 are the frequencies characterizing the anisotropic harmonic os- 
cillator potential. The energy levels for this system are simply 

111 

E n = (ni + -jtiuJi + (n 2 + -)hu} 2 + (n 3 + -)^3 , (16.2) 

where n = n 2 , n.3) with nj = 0, 1, (We have restored h here.) 

The spectrum is entirely discrete so that q = in the terminology of Sec. O. This 



means that BEC will not occur in the sense we have discussed so far; that is in the 
sense of symmetry breaking and as a phase transition. Nevertheless, as we will see, 
there is a temperature range which signals a point at which there is a sudden and 
dramatic increase in the number of particles in the ground state. If this criterion is 
adopted for BEC in place of the usual phase transition, then BEC will occur. The 
fact that a real physical system can exhibit such a behaviour is borne out by the 
recent experiments on cold alkali gases [44|-[47| . 



Because the atoms are neutral, we will use the particle number rather than charge 
as we have done earlier. The critical value of /1 is 

h 

H c = E = -{lo x +uj 2 + uj 3 ) . (16.3) 

There is an obvious similarity between this problem and BEC in a constant magnetic 
field, so that the methods used earlier should help us here. In order to present the 
analysis in as clear a way as possible we will first look at the case of the isotropic 
harmonic oscillator, returning to the anisotropic case later. 
Setting u)\ = u 2 = ^3 = oj in Eqs. |16.1| - |16.3| gives 



fie = \^ ■ (16.4) 

The g-potential defined in Sec. f| is 

00 00 00 

g-EEEM 1 - ze ~" En ] . ( 16 - 5 ) 

rii=0 ri2=0 rt3=0 

with z = e 13 ^ the fugacity. Because we are assuming an isotropic potential we may 

write the triple sum in Eq. pL6.5| as a single one. It is advantageous to introduce 



dimensionless variables x and e defined by 



x = P*u = ^,, (16-6) 
/x = • (16.7) 
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It is straightforward to show that 

oo 1 

q=-^2-{k + l){k + 2) In [1 - e ~ ik+e)x ] . (16.8) 



2 

fe=0 



There are several possible ways to proceed here. We will deal with the most direct 
way first. 

As we have done before, expand the logarithm in Eq. |16.8| in terms of its Taylor 
series to give 



* = E 2 {k + 1){k + 2) T, ■ ( 16 - 9 ) 

fc=0 n=l 

The sum over just involves a geometric series, and can be done with the result 

00 p — nex 

<7=E— (i-^r 3 • (le.io) 



n 

n=l 



The particle number is N = f3 1 (f^) • With T and uj fixed, we have from Eq. |16.7 
that 5 fx = —ftwSe. Thus 



N=- 1 (9q 



x 



(ff) =X^e" ne:c (l-e- nx )" 3 . (16.11) 



x n=l 

An equivalent way to obtain Eq. |16.11| is to begin with 



oo oo oo 

I -1 



ni=0 n2=0 rt3=0 

and using Eq. |16.2j (with all frequencies the same) follow through the same steps as 
led to Eq. |16.10| . The number of particles in the ground state is 

AT gr = ( e ex - I)" 1 . (16.13) 

Solving for e results in 

£= M i+ £>' (16 - 14) 

Since A^ < N we have 

e >iln(l + i). (16.15) 
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Figure 1: This shows e as a function of x — tkv/(kT). The solid curve shows the result found from 
a numerical evaluation of the exact harmonic oscillator sums for the isotropic case for N = 2000. 
The diamonds show the result of using our analytic approximation. 



This shows that e can only reach the critical value e = if x — > oo (corresponding to 
T — > 0), or if iV — > oo. For a finite particle number, as is the case in any experiment, 
e is always positive meaning that BEC as symmetry breaking will never occur. This 



is in agreement with the more general treatment in Sec. 

Even though e can never reach the value e = 0, if e gets very small then there 
can be a dramatic increase in the the number of particles in the ground state. In the 
actual experiments, u/(2ir) ~ 100 Hz, and T ~ fiK. In this case x defined in Eq. |16.6 
is small. The particle numbers of interest to the experiments are N ~ 10 3 — 10 6 . For 
a specific example we will pick iV = 2000 as in the rubidium experiment jH| and 
u/(27r) = 60 Hz. It is possible to solve Eq. |16. 11 numerically for e as a function of x 
for a given particle number. The result is shown in fig. 1. This calculation shows that 
over a very small range of x, corresponding to a very small range of temperature, e 
undergoes a very rapid decrease from large values to very small values. Associated 
with this sudden change in e is a sudden rise in the number of particles in the ground 
state. The ground state occupation number is shown in fig. 2. These numerical results 
show quite conclusively that BEC does occur in the sense that a significant number 
of particles occupy the ground state, and that the ground state occupation number 
increases very dramatically over a small temperature range. Because the behaviour 
of both e and N gT is smooth, there is no unique temperature which signals the onset 
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Figure 2: This shows the number of particles in the ground state as a function of x = ftio/(kT). 



of BEC. We will return to this point later. 

Rather than perform laborious numerical calculations we can try to obtain ap- 
proximate analytical results. Because the interesting behaviour is found for x « 1, 
the most obvious thing to try is to replace the sums over the discrete energy levels 
with integrals. However this is like regarding the energy levels as continuous rather 
than discrete. In the terminology of Sec. Ti this effectively replaces q = with q = 3, 
drastically altering the behaviour of the system unless sufficient care is taken in the 
conversion of sum to integral. We will return to this approach later. Instead of con- 
verting the sums to integrals, we will deal with the sums directly. For small x the sums 
in Eqs. |16.10| and |16.11| do not converge particularly rapidly, nor do they display in 
any transparent way the behaviour for small x. However we can use the Mellin-Barnes 
integral representation in Eq. |5.9| to great effect to obtain an asymptotic expansion 
valid for small x. 

Begin with Eq. |16.9| . It is possible to use Eq. [5]9] directly as in Ref. [48[ [|{|]; however 
as the ground state behaviour is of obvious interest it is convenient to separate it off. 
(It is just the k = term.) We will write 



q = q g r + q, 



16.16) 



66 



where 



£ 

n=l 



n 



-ln(l-e— ) , 



16.17) 



is the ground state contribution, and 



00 1 00 -n(fc+e)x 



k=l 



n=l 



n 



(16.18) 



represents the contribution from excited states. Using Eq. |5.9| on the RHS of Eq. |16.18| 
and performing the sums on n and k gives 



1 

2ri 



c+ioo 



daT(a)x~ a C R {a + 1) ( H (a - 2, 1 + e) 



+ (3 - 2e)( H (a - 1, 1 + e) + (1 - e)(2 - e)( H (a, 1 + e)J , (16.19) 

if we take c > 3 to avoid the poles of the integrand. If the contour is closed in the 
left hand plane, there are simple poles at a = 3, 2, 1, —1, —2, . . . and a double pole at 
a = 0. We find 



U4) 



.1 



(| " e)C*(3)aT 2 + i(l - c)(2 - e)( R (2)x- 1 + ■■■ . 



16.20) 



(The next term, which is of order lnx, is very complicated and will not be given here.) 
A similar procedure can be used to evaluate the particle number. Defining 



N = N gr + N ex , 



16.21) 



where iVg r was given in Eq. |16.13| and iV ex is the contribution from excited states, we 
find 



0?(3)*" 3 + 



e)CR®x- 2 - 



2x 



19 5 3 2 

h -e e 2 

12 2 2 



l-e)(2-e)(lnz + V(l + e)) + 0(1 
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It is possible to extend the calculations to any desired order in the expansion. 

We can now test how good the analytic approximation for the particle number 
is since we have already evaluated the exact expressions numerically. Let / be the 
fraction of the total number of particles in the ground state : 



fN. 



From Eq. |16.14| we have 



ex 



In 1 + 



fN) ' 



(16.23) 



16.24) 
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which determines e given x, f, N. If we keep only the first two terms in the expansion 



for N ex we have 



N^N gT + Ck(3)x~ 3 + (J - e)U2)af 



With Eqs. |16.23| and |16.24j this becomes 



(16.25) 



1_/)jVx 3 --Ch(2)x+ Ch(2)1b(1 + 



U3) 







16.26) 



The result of our approximate evaluation of e is shown as the diamonds in fig. 1. 
We have used the exact result for e to compute /, and given / and iV we can solve 
Eq. |16.26| for x. The approximate value for e is found using Eq. |16.24| . As can be 
seen from this figure the approximation is very good when e is small, but starts to 
break down when e becomes large. This is to be expected because by keeping only 
the first two terms in Eq. |16.22| we have treated e as small. (We have really assumed 
that e << x" 1 / 2 .) We will discuss another method for obtaining results with a wider 
range of applicability later. Even for e ~ 0.5 our approximate result is not too bad. 

We mentioned that a plausible approach if we want the behaviour for small x is 
to replace the sums over the discrete energy levels with integrals. This was done in 
Ref. [EH]. The result of this for the particle number is that only the first term in the 
expansion in Eq. |16.22| is obtained. This can be seen starting from 



~J> + !)(& + 2) 



-1 



16.27) 



k=l 



which follows from differentiation of Eq. |16.18| . Assume that there is a well defined 
transition temperature below which e = as for the normal unconfined Bose gas. 
From our discussion above, e really only vanishes for N oo which is often referred 
to as taking the thermodynamic limit. We prefer, following Ref. [28], to refer to 
this as the bulk approximation. The bulk approximation consists of taking e = in 
Eq. |16.27| , along with (k + l)(k + 2) « k 2 , and replacing the sum over k with an 
integral : 



N, 



bulk 



dkk 2 



-,kx 



l)- 1 = ( R (3)x- 



16.28) 



(A more systematic way of approximating the sum with an integral will be described 
later.) The bulk transition temperature is just the temperature at which N gr = 0. 
This gives 



a^buik 



" Cfl(3) ' 
N 



1/3 



16.29) 
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Fi gure 3: The ratio of the approximate to the exact ground state particle number is shown. The 
diamonds illustrate the result of using our approximation and the crosses denote the results found 
using the bulk approximation. 



Note that Xbuik = ^/(kT hu i k ), and that Eq. |16.29| holds for x > Xbuik corresponding 
to T < Tb u ik- The result of using the bulk approximation for the ground state particle 
number is that for x > Xbuik, 



N. 



bulk 



N -N. 



bulk 



N 



^bulk V 
X J 
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Fig. 3 shows a comparison of the bulk result with the approximation described above. 
The crosses show the ratio of the bulk approximation to the exact result for iV gr , 
and the diamonds show the ratio using N gT determined by Eq. |16.25| to the exact 
result. It is evident that the approximation in Eq. 16.25| is much better than the bulk 
approximation, especially as x becomes closer to the value where e undergoes its most 
rapid change. 

Having seen that it is possible to improve the bulk approximation for the ground 
state particle number significantly, we can see how Eq. |16.29| compares with the 
temperature quoted in the experiments. From Eq. 16.29| we have 

q l/3 

(16.31) 



bulk 



T 



N 



Ca(3) 



For rubidium with N = 2000 and Hlv/(2tt) = 60 Hz (which is the geometric mean of 
the frequencies) we find T bulk ~ 34.12nK, which is close to the value quoted in the 
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experiment [33]. Given the simplicity of the theoretical model this is perhaps a bit 
surprising. 

It is straightforward to see how the bulk temperature is altered by using the more 
accurate approximation. We know that Tb u ik is close to the experimental value, so 
assume that 



x = x bu ik(l +rj) 



(16.32) 



where rj is assumed small. Neglecting e in Eq. |16.25| we find 



Thus 



' 3 J + 2[C R (3)]2/3 
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bulk 
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(If we set / = in this result we obtain the result in Ref. |51], 52 1.) Ignoring /, 
which is typically a few percent, and taking iV = 2000 we find that the more accu- 
rate approximation reduces the temperature by around 10%. The larger the particle 
number the better the bulk approximation for the temperature becomes. However 
it should be clear from our study of the ground state particle number that the bulk 
approximation should be used with caution for other calculations. We will see that 
this is even more true when we study the specific heat. 

Before turning to the calculation of other thermodynamic expressions we will 
discuss a variation on the method presented above. This makes contact with the 
approach of ref. |k| which uses the Euler-Maclaurin summation formula to evalu- 
ate the exact sums defining the thermodynamic quantities. Previously we removed 
the ground state contribution from q. We will now remove the first excited state 
contribution as well. We may write 



<?ex = qi + q e 



16.35) 



where 



<?ex 



3 In 1 -e~ {l+e)x 



-(* + !)(* + 2) In 



k=2 



(16.36) 
(16.37) 



(See Eq. |16.8| .) We can expand the logarithm in Eq. |16.37| to obtain 



I ""^ g-n(fc+e)x 



k=2 



n=l 



n 



(16.38) 
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Now use the representation in Eq. 53, but only on part of the exponential. We may 
write 



-n{k+t)x 



g— n{h— l)a5g— n(l+e)as 



and use Eq. O on e""^- 1 )* only. This gives the contour integral representation 

1 



1 

2m 



c+ioo 

daY{a)x- a U l+a (e- {1+e)x 

c—ioo 

3 

+ oCr(« - 1) + Cfl(a 



16.39) 



where Li is the polylogarithm function introduced earlier in Eq. |6.21| . We take c > 3 
here to justify the interchange of sum and integral. The contour can be closed and 
the integral evaluated as usual. Lii +Q , is an analytic function of a, so that poles can 
come only from T(ot) or the Riemann (^-functions. We find 

9« = x- 3 Li 4 (e-( 1+ ^)+^^ 2 Li 3 ( e -( 1+ ^) 



+x- 1 Li 2 (e-( 1+ ^) +-ln 



-(l+e)x 



19 



240 



-x 



Jl+e)x _ 1 



+ ■ ■ 



(16.40) 



if the first few terms in the expansion are kept, and properties of the polylogarithm 
are used. A technical point which makes this approach useful for numerical com- 
putations is that as e becomes small the argument of the polylogarithms become 
approximately e~ x < 1. If we had not removed the first excited state then we would 
have obtained an expansion similar to that in Eq. |16.40| but where the argument of 
the polylogarithms was e~ ex ~ 1. We would then have found that terms in the ex- 
pansion in Eq. |16.40| would have become large as e got small, rendering the expansion 
inaccurate. Combining Eqs. |16.36| and |16.40| we find 

g cx = x- 3 Li 4 (e-( 1+ ^) +^- 2 Li 3 (e-( 1+ ^ 



+x- 1 U 2 (e- {1+e)x ) In 
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1 _ e -(l+e)s 
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This can be used as the basis for an accurate numerical approximation |53| , although 
some of the simplicity of the earlier purely analytic approximation is lost due to the 
presence of the polylogarithms which must be evaluated numerically. 
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17 Specific Heat for the Isotropic Harmonic Oscil- 
lator 



The specific heat has been used as a criterion for defining a critical temperature 
in finite systems [54]. A natural criterion is to define the critical temperature as 
the temperature at which the specific heat has its maximum. The behaviour of the 
specific heat for the unconfined Bose gas in three spatial dimensions is standard 
|I4] , |i~5| , |1~6| , [55]]. The specific heat is continuous with a maximum coinciding at the 
BEC transition temperature T c where the phase transition occurs. The derivative of 
the specific heat with respect to temperature is discontinuous at T c . The behaviour 

II. For D = 4 the 



of the specific heat in dimensions other than 3 is also known p7 



specific heat and its first derivative are continuous at the transition temperature T c , 
but the second derivative is discontinuous. For D > 5 the specific heat itself is not 
continuous at the transition temperature. In the cases of D — 1, 2 where no phase 
transition occurs the specific heat is a smooth function of temperature. We would 
expect that since no phase transition occurs for the Bose gas confined by the harmonic 
oscillator potential, the specific heat should be a smooth function of temperature. 

The starting point for an evaluation of the specific heat is the internal energy 
defined by 



U 



E4 



/3(E n -M) 



17.1) 



where the energy levels are given by Eq. |16.2j . For the isotropic harmonic oscillator 
we find 



Huj ^2 

k=0 



,(k+e)x 



-1 



17.2) 



when written in terms of the dimensionless variables x and e defined in Eqs. |16.6| and 
16. 7| . The term which involves 3/2 in Eq. |17.2| may be recognized as the total particle 



number. We have 

U _ 3 
ftw ~ 2 



N + - J2 k ( k + + 2 ) U (k+€)x - 1 



-i 



(17.3) 
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If we binomially expand the factor in square brackets in Eq. |17.3| we find 

jj oo oo 

— = i N + ^ E k ( k + x ) ( k + 2 ) E e ~ n(k+e)x 



_ 3 
JvIj ~ 2 
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The sum over k is easily performed with the result 

U 3 



fkjj 



-N + 3«i 



(17.5) 
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where 



00 _ 4 



ui = e~ n{1+t)x [I - e~ nx j . (17.6) 

n=l 



The specific heat is defined by 



Here u; plays the role of the volume for the free Bose gas. Because N is held fixed when 
we differentiate U, only the term 3u± will contribute to the specific heat. Because the 
temperature enters x = ftfrw, we have 

Care must be taken because N is fixed rather than e. From Eq. |16.11| we have 

N = e ~ neX ( l - e ~ nX ) > ( 17 - 9 ) 



n=l 



which defines e implicitly as a function of N and x. Differentiating both sides of 
Eq. |17.9| keeping N fixed gives us 

d(ex)\ S 2 
where 

S x = ^ne- nex (l -e~ nx ) , (17.11) 



71=1 
OO 



S 2 = ne- n{t+1)x (l - e" nx ) . (17.12) 

71=1 

Performing the differentiation in Eq. |17.8| and using Eq. |17.10| results in 

j = 3x 2 ^S 3 + S 2 -3^-j , (17.13) 
with 

00 j. 

s 3 = J2 ne ~ n(t+2)x (i - e ~ nx )~ ■ ( 17 - 14 ) 

71=1 
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Figure 4: The specific heat computed numerically for the isotropic harmonic oscillator is shown as 
the solid curve. The diamonds show the result using our approximation. The units for the specific 
heat are in factors of the Boltzmann constant k. The particle number is N = 2000. The maximum 
occurs for x ~ 0.0921. Our approximation breaks down for x below the point where the specific 
heat maximum occurs. 



This formula may be used to compute the specific heat numerically. The procedure 
used is to first solve Eq. |17.9| for e for a given N and x. Knowing e and x it is possible 
to evaluate the sums Si, 5*2 and S3 and hence determine the specific heat. The result 
of this calculation is shown in fig. 4 for N=2000. This figure shows that the specific 
heat is a perfectly smooth function of temperature in contrast to the result for the 
free Bose gas which has a discontinuous first derivative. A feature which the confined 
gas has in common with the free gas is that there is a well defined maximum. The 
maximum can be seen to occur at a value of x which is very close to the place where 
e undergoes its most rapid decrease and the population of the ground state starts to 
rise. We can therefore identify the specific heat maximum with the onset of BEC. 

Rather than perform detailed numerical calculations we can try to obtain an ana- 
lytical approximation for the specific heat by using the Mellin-Barnes representation 
to approximate the sums in Eq. |17.14| . If we assume x << 1 and e << 1 we find 

j ~ 12Ci?(4)x- 3 + % R (3)x- 2 + 2( R (2)x- 1 - l2e( R (3) X - 2 

-18e 2 ( R (2)( R (3)x~ 3 - 9e 2 a (3 ) x -4 + 9e 4 C£(3)ar 6 . (17.15) 



The result of using Eq. |17.15| to calculate the specific heat is shown as the diamonds 
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in fig. 4. It can be seen that at temperatures below the point where the specific heat 
has its maximum the approximation is quite good. For values of x < x m ~ 0.0921 
the assumption that e is small is no longer valid and the simple approximation breaks 
down. It is possible to obtain a better approximation using polylogarithms. 

We can also see what happens if the bulk approximation is used. From Eq. |17.4 
if we replace the sum on k with an integral and use k(k + l)(k + 2) « k 3 we obtain 



bulk 



1 r 00 
N + - e ' neX / dkk 

n=l ^° 



3g _ nkx 



-iy + 3x- 4 Li 4 (e- 



For T < Tb u ik we have e = so 



gives the internal energy. The specific heat for T < Tbuik is 

= 12k( R (4)x- 3 . 



°bulk 



dT 



17.16) 



(17.17) 



17.18) 



N,LU 



This is observed to correspond to only the first term on the right hand side of 

Eq. irrrq . 

The bulk approximation for the specific heat when T > Tb u ik is more involved. 
We first need the number of particles when e ^ 0. This can be obtained exactly as in 
Eq. |16.27| by approximating the sum over k with an integral. The result is 



N hulk = x~ 3 Li 3 (e— ) . 
We define Tbuik as before by setting e = in Eq. |1719| and obtaining 



Eliminating iV between Eq. 17.1S and Eq. |17.20| results in 



Li 3 (e" ex ) = Ca(3) 
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bulk 



We can now use Eq. |17.16| to give 



^buik — 3 A; 



k 



d_ 
df 



{T 4 Li 4 (e— )} 



By differentiating the polylogarithm it is easy to see that 



d(ex) 
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Li 3 (e" 
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Differentiating Eq. 17.21| with respect to T holding N and uo fixed results in 

'd{ex)\ _ 3Ci ? (3)T b 3 ulk T- 4 



dT 



[17.24) 



Putting these results together yields 
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For comparison, we can use Eq. |17.20| in Eq. |17.18| to obtain 
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fciV "C«(3) V^buik, 
^buik in Eqs. |17.25| and |17.26| we find 
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The bulk approximation leads to a discontinuous result for the specific heat. This 
contrasts with both the free Bose gas, which had a continuous but not smooth be- 
haviour, and the exact result which was shown in fig. 4 which is continuous and 
smooth. Although Xb u ik is close to the BEC temperature determined by the true 
specific heat maximum, the bulk approximation for the specific heat shows features 
not present in the true system. 

Leaving aside the fact that the bulk approximation for the specific heat is dis- 
continuous, we can see how good an approximation it is when compared with the 
true result. This is shown in fig. 5. The crosses show Cb u ik/Ccxact and the diamonds 
show the ratio of the approximation in Eq. |17.15| to the exact result. Although the 
agreement between the bulk value and the exact value is good near the specific heat 
maximum, it is off by about 25% once x has reached 0.4. In contrast Eq. 17.15 is 
within about 1% of the true value at this point. 



18 Anisotropic Harmonic Oscillator Potential 

The case where all of the oscillator frequencies are different is more difficult to analyze 
than the isotropic case. Nevertheless the Mellin-Barnes representation can still be 
used to great effect H4~9"|. We have 



E n = %(n\Ux + n 2 u)2 + ^3^3) + ii c 
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Figure 5: The diamonds show the ratio of our approximation for the specific heat to the exact 
result. The crosses denote the ratio of the bulk specific heat to the exact value. N — 2000 is taken. 
Both results become increasingly inaccurate below the specific heat maximum. 



if Eqs. |16.2| and |16.3| are used. We will define dimensionless variables x\, x%, X3 by 



analogously to x in Eq. |16.6| . We will also define 

to be the average frequency, and 

a = h(3Vt = ]-(x 1 + x 2 + x 3 ) 

Finally we will define 
analogously to Eq. \16.7[ We then have 



(18.2) 



18.3) 



18.4) 



18.5) 



f3(E n — /i) = riiXi + ae . 



i=l 



(18-6) 
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The g-potential in Eq. |16.5| reads 



oo oo oo 



g = -^^^lll[l- e -{n 1 x 1 +n 2 x 2 +n 3X3+a e)^ ? ^g ^ 

ni=0 n 2 =0 n 3 =0 

when it is written in terms of dimensionless variables defined above. 
As before, we expand the logarithm in its Taylor series 

oo ., oo oo oo 

e -n(n 1 x 1 +n 2 X2+n 3 x 3 +ae) Qg g^ 



^*_< ^ 

71=1 ?ii=0 ri2=0 ri3=0 
oo .. 

£ -e- nae (l - e^ 1 ) -1 ^ - e-™ 2 )- 1 ^ - e"^ 3 )" 1 , (18.9) 



n 



where the second line has followed from performing the sums over n 1; n 2l n 3 . We can 
separate off the ground state from the sums as we did in the isotropic case, but this 
gets slightly messy here so is left as an exercise for the reader. Instead we will apply 
Eq. |5.9| directly to the exponential in Eq. |18.8| to obtain 

g=— / dzT(z)C R (z+l)C B (z,ae\ X ) , (18.10) 

where (b is the three dimensional Barnes (^-function defined in Eq. [11 .25| . This func- 
tion is analytic for 3l(z) > 3 with simple poles at z = 3, 2, 1. Using the values of the 
residues at these poles given in Eqs. 11.28| - T1.30 we find 

q ~ + --e )a 

0l ' (2) [Qa 2 e 2 -18a 2 e + 9a 2 + x l x 2 + x 2 x 3 + x 3 x 1 ] + --- . (18.11) 



12 X1X2X3 
The particle number is 



N = r l (P) =—(%') , (18.12) 



if we use the fact that fixed T and u>t corresponds to fixed Xj, and hence fixed Q and 
a; thus, <5/i = — hQ5e = —i3~ 1 a5e. From Eq. |18.8j we find easily that 
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C(3) + e Wl + .... (18.15) 
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The first two terms on the right hand side of Eq. 18.15| are easy to obtain, and in fact 



follow directly from differentiating the expression for q in Eq. |18.11| . The third term 



(ae) -1 requires some explanation. Its origin is the pole of the integrand of Eq. |18.14j 
at z — 1 which is a double pole since (r(z) — l/(z — 1) + 7 + • ■ ■ near z — 1 and 
has a pole at z — 1. To evaluate the third term completely we need to know the finite 
part of (b at z = 1 in addition to the pole part. This is not known to my knowledge; 
however, if we are interested simply in the case where ea is small, which is the case 
on physical grounds, then we may use the fact that (b{z, a|x) ~ a~ z gives the leading 
contribution (coming from the n = term in the sum). Keeping only the dominant 
term for small ae results in Eq. pL8. 15| . Another way to check the validity of Eq. |18. 15 



is to remove the ground and first excited state from the sums in Eq. |18.13| in the same 
way as we did for the isotropic oscillator, and express N as in Eq. |18.14| but where 
polylogarithms arise in place of the Riemann (^-function. This avoids the problem of 
the double pole at z — 1 because the polylogarithm is analytic there. Furthermore, 
this allows the limit e — > to be taken because the argument of the polylogarithm 
involves e~ Xmin where x m j n is the smallest of the Xj. 

A calculation similar to that outlined for the isotropic harmonic oscillator may be 



performed here. The details of this may be found in Ref. ||49|| . Again the results of 
the analytic approximation may be compared with the exact results. The agreement 
between the two is as good as in the isotropic case. For N = 2000 and uj\ = 0J2 = 
2407r/v / 8 s -1 and U3 = 2A0n s -1 the specific heat maximum is found numerically to 
occur at a = 0.106. Using Eq. |18.4j this may be shown to correspond to a temperature 
of 30.9 nK which can be compared to 31.3 nK if we approximate the system by an 
isotropic harmonic oscillator whose frequency is the geometric mean of the three 
frequencies, and to 34.1 nK if the bulk approximation is used. The fact that the 
oscillator is anisotropic only changes the temperature by a few percent from the 
simpler isotropic calculation. 

A final point we wish to discuss is that the (ae) -1 term in Eq. 18.15 is observed 
to be the first term in the expansion of 



jV gr = ( e M - l)" 1 , (18.16) 



in powers of ae. Thus we can interpret the first two terms in Eq. |18.15| as the 
contributions from the excited states : 

». S M + M(' ( )„. 08.17) 



XiX 2 X 3 XiX 2 X 3 V2 

An analysis similar to that presented for the isotropic harmonic oscillator may now 
be given. 

N cx ~ I*®- + , (18.18) 
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if we assume that e << 3/2. By saying that at the onset of BEC we have N e - 
we obtain the equation 



N 



N 



Cr{3) 3Cr(2) 

+ ot 



;i8.i9) 



which determines T. (Setting a = in Eq. |18.19| is the bulk approximation.) Using 
Eqs. |18.2| and |18.3| results in the cubic equation 



where 



N ^aT 3 + bT l 



^ 3 C«(3) 



2fb 2 UJiUJ2i02, 

Since N is large and T is small we may assume 



(18.20) 

(18.21) 
(18.22) 
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with 77 small. It is straightforward to obtain 



kT « huj 
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^(2) 
2uCf (3) 



(18.24) 



where 



,1/3 



As before, the second term in Eq. |18.24j represents the 



correction to the bulk result. Because we have assumed N ex = N, there will be a 
small correction to this if we compare the result with the specific heat maximum due 
to the fact that at the specific heat maximum there is a non-zero occupation for the 
ground state particle number. Taking iV = 2000 and the frequencies quoted before 
gives T ~ 31.9 nK which is about a 6% reduction from the bulk value. More details 



on the anisotropic oscillator may be found in Ref. 49 . 



19 Density of States Method 

We have already discussed the bulk approximation in which sums are replaced with 
integrals. As a rough approximation for calculating the BEC transition tempera- 
ture this might be sufficient; however we saw that there was a more accurate way of 
obtaining approximate results. In Refs. |)T, Q it was discussed how an improved 
approximation in which sums were converted to integrals could be achieved by includ- 
ing a modified density of states factor for the harmonic oscillator. For the isotropic 



80 



oscillator we had the energy levels E k = (k + 3/2)fkv with degeneracy \(k + 1) (fc + 2). 



(See Sec. [L6|.) Grossmann and Holthaus suggested that the sum over k could be con- 
verted to an integral over k but with a density of states factor |fc 2 + |/c in place of 
simply jk 2 as in the bulk approximation. By a change of variable from k to E, the 
density of states factor becomes 

/ „n 1 E 2 3 E , 

For the anisotropic oscillator the situation is more complicated, and the density of 
states was parametrized by 

^ = (19 ' 2) 



in Refs. pi], |52|j where 7 is a constant which depends on the frequencies and which 



had to be evaluated numerically. (In Eq. |19.1 uo is the geometric mean of the three 



frequencies. It was later shown [^g, |57j how 7 could be arrived at analytically, and it 
is this approach we will discuss here. Furthermore it is possible to set the problem of 
conversion of the sums to integrals in a much wider context leading to other possible 
applications. 

Let H be the Hamiltonian operator for some system with E n the energy levels. 
Let M{E) be the number of energy levels for which E n < E. We can write 

N{E) = Y J G{E-E n ) , (19.3) 

n 

where 8(x) is the Heaviside distribution (or step function) defined by 

( 1, x > 

9(x) = { 0, x < . (19.4) 




The aim is to treat the energy levels, or at least some part of the energy spectrum, 
as a continuous set rather than a discrete set by introducing the density of states 

p(E)dE = Af(E + dE)-Af(E) 

Since 8'(x) = 5(x) we can write 

p(E) = J2^(E-E n ) . (19.6) 



Now define 



POO 

K{t) = / dEe- tE p(E) , (19.7) 
Jo 
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which is just the Laplace transform of the density of states. If we use Eq. 19.6 , then 
K(t) becomes 



Kit) = e~ tEn , 

n 

The inversion formula for Laplace transforms gives 

i pc+ioo 

P(E) = — dte tE K(t) 
2m J c _ ioo 



:i9.8) 



(19.9) 



for c G M with c > 0. 

The importance of Eq. |19.7| or Eq. |19.9| is that if we can evaluate K(t) given in 



Eq. p.9.8| then the density of states may be found. The problem of evaluating Eq. |19.8 



for a given E n is sometimes not too difficult, as in the case of the harmonic oscillator. 
We will do this in Sec. |19.1| below. 

Even if it is not possible to perform the sum in Eq. |19.8| exactly, it is still possible 
to obtain a knowledge of the density of states. This is because the density of states is 
determined by the asymptotic behaviour of K(t) as t — > which is known for a wide 
class of operators. We will assume that as t —>■ 



K(t) ~ J2 Cit ' n + °( r 



19.10) 



8=1 



for some coefficients q and powers with r\ > r 2 > 

1 



r k > 0. Noting that 



r 



Tin 



(19.11) 



for Ti > it is easy to see from Eq. |19.7| that 

k 



n-1 



i=l 



(19.12) 



Thus a knowledge of the coefficients q in Eq. |19.10| gives us the density of states. 
(The result may be extended to < but it requires a more powerful method than 
Laplace transforms [K? 



19.1 Application to the Simple Harmonic Oscillator 

We will show how the density of states method may be used to obtain results which 
agree with those found in Sec. [16]. The g-potential is 

q = q gT + g ex , (19.13) 



82 



where 

g gr = - In (1 - ze- pEo ) , (19.14) 
is the lowest mode or ground state contribution, and 

q ex = -J2 ln i 1 - ze ~ 0En ) > ( 19 - 15 ) 

n>0 

represents the contribution from the excited states. We can write 

q cx = £ £ e"^"^ , (19.16) 

k=l ' n>0 

since fi c = E . We now convert the sum over n in Eq. |19.16| into an integral with the 
use of the density of states : 

~ E / dEp{E)e~ kfiE • (19.17) 

fc=l fc -^-Bo 

The lower limit on the integral corresponds to the first term in the sum over n in 
Eq. [T9~T6 . 



We now need p(E). To obtain this we first calculate K(t) defined by 

K(t) =J2 e ~ t(En ~ Eo) ■ ( 19 - 18 ) 



n>0 



Because we only want the behaviour of K(t) for negative powers of t we can add back 
the n = term to obtain 



oo oo oo 



K(t) = ^ ^ ^ ^ ^ ^ e -th{n 1 ui 1 +n2LU2+n 3 u) a ) _ j 

ni=0 ri2=0 ri3=0 
3 

= JJ (1 - e - *^)" 1 - 1 . (19.19) 

3=1 

It is now straightforward to expand K(t) in powers of t using 

(l-O-i = i + i + |L + 0(x 3 ). (19.20) 



(This expansion just involves the Bernoulli numbers |I0"|| .) We find 

K{t) ~ Cl r 3 + c 2 r 2 + cgt- 1 + 0(1) , (19.21) 
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where 



c 3 



1 



1 1 1 
+ + 



2ft? \ u)\U)i UJ2UJ3 UJ3UJ1 



1 



^3 ^2 ^1 



12/1 \a;ia;2 ^1^3 0*2^3 
We can now use Eq. |19.12| to find 



3 



3 

^2 



3 

u; 3 



(19.22) 
(19.23) 

(19.24) 



p(£) 



-C1E 2 + c 2 E + c 3 . 



(19.25) 



The first two terms agree with the ansatz of Ref. |5T], |52] and actually represent a 
direct evaluation of their coefficient 7. It is easy to see by comparing Eq. |19.2| with 
Eq. |ilT25| that 



(19.26) 



KM" 



7 = -{UJ X + U 2 + ^ 3 )(^i^ 2 ^3) 1/3 • 

Finally we note that for the isotropic oscillator we have c\ = (ftu)~ 3 ,C2 
and c 3 = (ftu)^ 1 which agrees with a direct conversion of the degeneracy factor 
\(k+l)(k + 2). 

We now return to the evaluation of q cx in Eq. |19.17| using Eq. |19.25| . Before doing 
the integrals we will make the simplifying assumption that j3(Ei — Eq) << 1 so that 
we can replace the lower limit on the integral with zero. We then have 



E 

k=l 

00 

E 

k=l 



k=l 

00 k^Kle 



J dE^ciE 2 +c 2 E + c 3 y- 



k/3E 



k 



Cl 



+ 



c 2 



+ 



c 3 



(k(3) 3 (k(3) 2 (kp) 



(19.27) 



if we use Eq. |18.5| with \i c = E = hh(u>i + u 2 + cu 3 ) = |M1 The sums in Eq. |19.27 
are just polylogarithms, and we obtain 



Cl /5- 3 Li 4 (e-^) + c 2 /r 2 Li 3 (e-^) 
+c 3 /5- 1 Li 2 (e- /3 ^) . 



(19.28) 



The expansion of the RHS in powers of a = (3ftQ generates Eq. |18.11| . The resulting 
discussion is the same as before. 



19.2 Application to BEC in Cavities 

Another problem which can be studied is the situation where a system of spin-0 
bosons is confined to a cavity of arbitrary shape. The results for the g-potential in 
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Eqs. |19.13| - |H).17| still apply, but of course the density of states factor will be different 
than for the harmonic oscillator potential. We can concentrate on the kernel function 
defined in Eq. |19.25| as before and evaluate the coefficients q appearing in Eq. |19.10| 
as t — ► 0. General spectral analysis |0| gives 



K(t) ~ (47rr)- D / 2 |a + ^r 1 / 2 + O(r)} , (19.29) 

as r — > 0, where we have defined 

K(r) = tr ^ e -(-v 2 +v(x)) j (lg 3Q) 

The coefficients ao and a% appearing in the asymptotic expansion Eq. pL9.29| are 

a = Vv, (19.31) 
ai = ^dVv, (19.32) 

where <9Vs denotes the volume of the boundary of the cavity, and b = — 1 for Dirichlet 
(ie fields vanish on <9Vx) and b = +1 for Neumann (ie normal derivative of fields vanish 
on dVs). The boundary of the cavity has been assumed to be smooth and we take 
the spatial dimension D > 2. 

The case of interest to us corresponds to r = ^- so that the operator appearing 



in Eq. |19.30| is the Hamiltonian. The coefficients appearing in Eq. p.9.10 



arc 



ci = ^(^J (with n = D/2) (19.33) 

b / m \ P-l)/2 
c 2 = fV s (—) (with n — (D— l)/2) . (19.34) 

We therefore obtain from Eq. |19.12| the density of states 

K FfTJ/Sf D/2 " + i#) £< " 3,/2 ' (19 - 35) 

to the order we are working. The first term on the right hand side is sometimes called 
Weyl's theorem. (See Ref. |5"9" , |BT| for example.) From Eq. |19.17| , again approximating 
the lower limit on the integral with 0, gives 

qc , - Y,-^—{ c ^y DI2 + c ^y {D ' l) ' 2 } 



k=l 



V^\ T D U D/2+l {e^^) + \dVz\ T {D - l) U {D+l)/2 {e^-^) , (19.36) 
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where 



A T = V ^, (19.37) 
V m 

defines the thermal wavelength. From this it can be seen that the expansion parameter 
is Xt/L where L = is a typical length scale associated with the cavity. In the 

infinite volume limit we can ignore the term due to the boundary and the standard 
result for the free Bose gas found in Sec. ^T] is recovered. 

We are now in a position to see how finite volume effects can modify the tem- 
perature associated with BEC in a finite cavity. For a system confined to a finite 
volume the energy levels will be discrete in general, so that the criterion for BEC 



as symmetry breaking discussed in Sec. [14] is not met and symmetry breaking and a 



first order phase transition will not occur. However in the bulk limit (ie the infinite 
volume limit) if D > 3 we do have BEC as symmetry breaking. Therefore on physical 
grounds we expect that for a very large, but finite, volume the specific heat should be 
smooth, but should approach the non-smooth behaviour as the infinite volume limit 
is taken. For large volumes, although the chemical potential will not be able to reach 
its critical value, it can become arbitrarily close much as in the harmonic oscillator 



potential discussed in Sec. 16. This should correspond to a sudden rise in the number 



of particles in the ground state 
a_ 
dfi 



Using N cx = (3 1 jf-q cx with q cx given in Eq. |19.36| results in 



N ex ~ V s A^Li D/2 (e^-^) + ^dV s \- {D - 1] U {D „ 1)/2 (e^~^) . (19.38) 
The ground state particle number is 

jV gr = [ePfa-d - l]" 1 , (19.39) 
as usual. Our expansion parameter is £ = L/Xt with £ >> 1. Let 

dV^ = kL - 1 , (19.40) 

where k is a dimensionless constant which depends only on the shape of the boundary. 
If we let / be the fraction of particles in the ground state at the specific heat maximum, 
then 

/5(/i c -/i)=ln(l + -^) . (19.41) 
Assuming that fN » 1 we have 

/?(^ - ~ -L « 1 . (19.42) 
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We will approximate 

l t)K l 

N ex ~ Z D Li D/2 (e-7N) + _^-iLi (D _ 1)/2 (e-7iv) . (19.43) 

For .D > 3 the order of the polylogarithm is larger than 1 and we may further simplify 
this expression to 

N cx * eU§) + ^^U^Y 1 ) ■ (19.44) 
In the case of D = 3 we use Lii(^) = — ln(l — z) and approximate 

Lii(e~7^) ~ ln(/iV) . (19.45) 

Thus for D = 3 we have 

Ne^eU^ + ^eMfN). (19.46) 

In either case we have 

N cx = N-N gI = (l-f)N . (19.47) 
For D = 3 Eq. p^46| and Eq. ggg^ give us 



(l-f)N^ ( R (^)e + b -^e HfN) , (19.48) 
resulting in a cubic equation for £. The bulk value of £, call it £ , is given by 

N = Cr^o ■ (19-49) 

We are assuming here that the boundary corrections to the bulk result are small, so 
that £ in Eq. |19.48j should be close to the bulk value £ . Write 

£ = £o(l + 77), (19.50) 
where 77 is small. Working to first order in small quantities we obtain 

(1 " f)N * U§)$(1 + 377) + ^£ 2 HfN) . (19.51) 



It is easily seen that 



± (19 52) 



If T denotes the bulk temperature, then 



T ~ T (l + 277) . 



19.53) 



(This follows trivially from Eq. |19.50| since £ oc A^ 1 oc T 1 / 2 .) In the case of Neumann 
boundary conditions where b = +1 we can conclude that the temperature is lowered 
from the bulk value. For Dirichlet boundary conditions the boundary correction to 
the bulk result tends to increase the temperature, but the overall effect depends upon 
the relative size of the two terms in Eq. |19.52| . 

For spatial dimensions D > 3 we have in place of Eq. |19.51 



(l-f)N 
The bulk value of £ is given by 



bK D 

T C*(— 



D-l 



N = C 



R 



Writing £ as in Eq. |19.50| and working to first order in rj gives 



V = — FT 



f_ 

D 



nb ( t 



D-l' 



-l/D 



N' 



-l/D 



19.54) 



19.55) 



(19.56) 



Eq. |19.53| still relates T to the bulk temperature. 

BEC in a rectangular cavity has received numerous treatments. (See Refs. [54 
for a review of the earlier work, and Ref. [Ml for a recent analysis.) More details on 



the approach presented here can be found in Ref. [57 
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